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ABSTRACT 


A  study  was  performed  to  evaluate*  the  effect  of  concurrent  photon 
amplification  (CPA)  on  Kodak  Tri-X  Pan  emulsion  in  regard  to  image 
quality.  The  study  evaluated  detective  quantum  efficiency  (DQE), 
modulation  transfer  function  (MTF ) ,  and  information  content  (IC)  as  a 
function  of  exposure  in  a  comparison  between  normal  exposure  and  CPA. 

The  results  showed  CPA  to  enhance  image  quality  relative  to  normal 
exposure  at  very  low  exposure  levels.  As  the  exposure  was  Increased, 
normally  exposeG  images  became  comparable  and  then  exceeded  CPA 
exposures  in  image  quality.  The  exposure  level  at  which  normal  exposure 
became  preferential  was  distinct  with  each  image  quality  measure. 
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INTRODUCTION 


The  purpose  of  this  section  is  to  give  a  background  to  the 
experiment.  The  historic  development  and  mechanism  of  concurrent  photon 
amplification  (CPA)  will  be  discussed.  Measures  of  image  quality  will 
be  defined  and  the  effect  of  CPA  on  each  measure  will  be  postulated. 
Finally,  some  measurement  techniques  of  image  quality  appropriate  to 
the  laboratory  will  be  derived  or  defined. 

HISTORY 

The  development  of  CPA  can  not  be  attributed  solely  to  one 

individual  or  group.  In  the  lS30's,  it  was  known  that  some  additional 

light  sensitivity  in  silver  halide  films  could  be  gained  by  a  slight 

fogging  exposure  near  the  time  of  the  imaging  exposure^.  Fogging  was 

achieved  by  very  briefly  exposing  the  film  to  ambient  light,  then 

placing  the  film  into  the  camera  and  making  the  imaging  exposure.  The 

pioneer  of  making  the  fogging  exposure  coincident  with  the  imaging 

2  4 

exposure  was  Cole  in  1S72  ’  .  In  contrast  to  the  primitive  technique 
of  opening  the  film  back  to  ambient  light.  Cole  designed  and  built  a 
device  using  compact  electronics  and  light  emitting  diodes.  His  CPA 
device  would  synchronize  a  controlled  exposure  from  the  diodes  with  the 
operation  of  the  shutter  producing  a  predictable  level  of  fogging  and 
light  sensitivity  enhancement.  His  device  was  successfully  marketed  as 
a  modification  to  a  number  of  existing  cameras  and  enjoyed  some  success 
among  individuals  interested  in  low  light  photography.  To  further 
appreciate  the  success  of  his  product,  an  understanding  of  the  CPA 
mechanism  and  its  effect  is  needed. 
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MECHANISM 

The  light  sensitive  portion  of  a  silver  halide  emulsion  consists  of 
a  solution  of  silver  halide  crystals  suspended  in  a  gelatin  emulsion. 
One  photon  of  sufficient  energy  incident  upon  the  crystal  will  elevate 
one  Ag+  ion  to  an  Ag°.  The  Ag°  is  unstatle  as  a  singular  entity 
and  will  normally  revert  to  the  Ag+  state.  However,  if  four  or  more 
Ag°  appear  as  a  localized  aggregation,  then  all  four  will  be  stable 
as  metallic  silver.  The  four  Agc  aggregate  site  is  known  as  a  latent 

image  site  and  is  capable  of  being  a  catalyst  for  the  chemical 

reduction  of  the  remainder  of  the  crystal  during  development.  The 
formation  anc  stability  of  the  latent  image  site  is  known  as  the 

Gurney-Mott  mechanism1^.  The  purpose  of  CPA  is  to  provide  a  portion 

of  those  four  photons  generally  across  the  emulsion.  Anywhere  the 
remaining  necessary  photons  from  the  imaging  exposure  strike  a  crystal, 
a  latent  image  will  form. 


FILM  SPEED 

Because  fewer  photons  are  needed  during  the  imaging  exposure,  there 

is  an  apparent  film  speed  increase.  Speed  increases  of  six  stops  have 

4 

been  reported  using  Tri-X  Pan  emulsion  .  This  represents  a  very 
significant  enhancement  but  begs  the  question:  What  effect  does  CPA 
have  cn  image  quality?  One  answer  is  obvious,  if  there  is  an  image 
where  there  was  none  before,  there  is  an  improvement  in  image  quality. 
To  evaluate  that  improvement,  some  measures  of  image  quality  must  he 
defined  or  derivec. 
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CONTRAST 

The  effect  of  CPA  on  contrast  and  the  value  of  contrast  as  an  image 

quality  indicator  are  both  very  significant.  The  change  in  density  vs 

change  in  log  exposure  is  called  the  contrast  or  gamma  in  a 

photographic  system.  A  gamma  that  is  too  low  to  provide  an  acceptable 

presentation  of  the  scene  content  will  render  an  unacceptable  image. 

This  is  due  to  a  certain  minimum  level  of  change  in  density  that  can  be 

visually  or  instrumentally  detected.  A  gamma  equal  to  1.00  is  commonly 
3 

desired.  Due  to  the  nature  of  film,  the  gamma  for  the  mid-region  of 
exposure  will  be  higher  than  the  gamma  in  either  the  tee  (extreme  low) 
or  shoulder  (extreme  high)  region  of  exposure.  The  negative  deviation 
of  gamma  from  this  desired  level  degrades  the  image  quality.  CPA  will 
tend  to  increase  the  minimum  film  density,  thereby  reducing  the  gamma. 
This  effect  will  continue  until  the  irradiance  from  the  CPA  device 
becomes  insignificant  compared  to  the  scene  irradiance.  At  that  time 
the  gamma  will  increase  to  the  gamma  normally  held  by  the  film. 

NOISE 

Noise  in  a  photographic  system  is  the  disturbance  of  the  image 
caused  by  the  granular  nature  of  the  silver  halide  emulsion.  A  standard 

measure  of  emulsion  grain  noise  is  granularity,  which  is  proportional 

—  2 

to  the  root  mean  square  change  in  density,  d(D)  or  density  variance 

about  some  mean  level.  This  is  usually  measured  with  a  small  aperture 

densitometer  at  a  given  overall  density  level.  The  industry  standard, 

set  by  Kodak,  is  to  measure  this  density  variance  at  an  overall  density 

5 

of  1.00  with  a  48  micron  circular  aperture  scanning  densitometer. 
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If  the  number  of  photons  falling  upon  the  emulsion  follow  a  Poisson 
distribution,  then  the  distribution  of  the  developed  grains  should  also 
follow  the  same  form.  In  a  Poisson  distribution,  the  variance  will 
increase  with  the  increasing  number  of  developed  grains.  Density  is 
defined  as, 

density  =  -  log( transmission) ,  (1) 

where 

transmission  =  1  -exp,  (2) 

and  exp  =  the  fraction  of  grains  that  are  exposed.  For  the  domain  where 
the  Poisson  variance  equals  the  number  exposed  and  the  number  exposed 
is  linearly  increasing  with  increasing  exposure,  density  variance  will 
increase  with  increasing  exposure. 

The  impact  of  CPA  upon  the  ncise  level  in  the  emulsion  comes 
through  the  generation  of  fog.  By  its  nature,  CPA  will  generate  an 
increased  level  of  fog  relative  to  normal  exposures  and  by  association 
ar,  increase  in  the  noise  level  of  the  emulsion.  This  is  related  to  the 
quantum  sensitivity  distribution  of  the  emulsion  which  is  the  relative 
number  of  grains  that  become  developable.  Crystals  which  require  four 
photons  to  achieve  a  latent  image  site  are  said  to  have  a  quantum 
sensitivity  of  four.  In  this  situation,  CPA  could  provide  li%  of  the 
required  irradiance  to  create  the  latent  image.  For  crystals  *  ith  a 
much  higher  quantum  sensitivity,  e.g.,  20,  CPA  could  provide  a  higher 
percentage  of  the  photcr.s  necessary.  Measured  sensitivity  distributions 
of  various  workers^  using  special  test  emulsions  range  from  a  minimum 
level  of  3  to  over  20,  and  are  shown  in  figure  1. 
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figure  1,  Measured  Sensitivity  Distributions 
These  results  are  not  contrary  to  the  classic  Gurney-Mott  mechanism 
because  these  are  considered  probablistic  processes. 

The  sensitivity  distribution  of  the  emulsion  shows  both  the 
potential  benefit  and  problems  of  CPA.  A  level  of  CPA  that  would  give 
optimum  enhancement  for  grains  requiring  high  quantum  exposure  would 
generate  fog  in  grains  needing  low  quantum  exposure.  The  compromise  is 
immediately  apparent,  the  more  the  image  is  enhanced  with  CPA,  the  more 
fog  and  noise  generated. 


MODULATION 

Modulation  is  defined  by  the  equation 

modulation  =  (Dmax-Dmin)/(Dmax+Dmin).  (3) 

Where  Dmax  is  the  maximum  density  in  the  region  of  the  boundary,  say 
the  dark  side  of  an  edge.  Dmin  is  the  minimum  density  on  the  alternate 
side  of  the  boundary.  CPA  will  raise  the  mean  level  of  both  Dmin  and 
Dmax  in  a  non-linear  manner,  therefore  the  overall  effect  is  uncertain. 
When  the  change  in  Dmax  is  greater  than  the  change  in  Dmin,  there  will 
be  an  improvement  in  modulation.  However,  Dmax  will  reach  some 
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ultimately  limiting  value  and  while  Dmin  continues  to  rise,  CPA  will 
have  the  effect  of  reducing  modulation.  If  the  spatial  relationship 

between  Dmax  and  Dmin  is  varied,  and  modulation  measured  at  each 

spatial  frequency,  the  resulting  plot  of  modulations  vs  spatial 

frequency  constitutes  the  familiar  modulation  transfer  function  (KTF). 
KTF  is  a  very  valuable  tool  to  evaluate  the  information  carrying 
capability  of  an  emulsion.  A  specific  means  of  determining  KTF  will  be 
discussed  later. 

The  limit  of  resolution  is  the  spatial  frequency  where  the  minimum 
discernable  modulation  occurs.  This  minimum  discernable  modulation  may 
be  determined  by  human  or  instrumental  means.  For  both  CPA  and  normal 
exposure  it  is  possible  to  predict  the  effect  of  increasing  exposure 
upon  resolution  in  the  exposure  region  when  the  change  in  Dmax  is 

greater  that  the  change  in  Dmin.  A  relationship  can  be  derived  between 
minimum  discernable  modulation  anc  resolution  in  terms  of  exposure7. 
Their  relation  can  be  seen  by  considering  two  adjacent  squares 
receiving  EX1  and  EX2  levels  of  exposure  in  terms  of  photons  per  area. 
These  two  squares  are  discernable  when 

EX1  +  k(EXl)1/2=  EX2  -  k(EX2)1/2,  (4) 
or 

E XI  -  EX2  =  k  (  (EX1)1/2  +  (EX2)1/2).  (5) 

The  value  of  k  used  is  a  function  of  the  confidence  required  to  insure 

a  discernable  difference  exists.  Altman  and  Zweig  have  shown  that  if  k 

21 

=  5,  there  is  a  1  part  per  million  error  rate.  Therefore,  k  is 
normally  set  equal  tc  5.  A  higher  tolerable  error  rate  would  result  in 
a  lower  value  of  k.  These  two  squares  are  considered  as  part  of  a 
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larger  image  with  an  overall  average  exposure  of  EXB.  Modulation  is 


defined  as 


M  =  (EX2  -  EX1)/(EX2  +  EX1),  (6) 


EX1  4  EX2  =  2EXE  (7) 


(EX1)1/2  4  (EX2)1/2  =  2  (EXB)1/2.  (8) 

The  combination  of  equations  E,7,  and  8  into  6  leads  to  an  expression 
for  minimum  discernable  modulation, 

Kmin  =  k(EXB)"1/2.  (9) 

If  each  patch  has  an  equal  area, . denoted  area,  the  resolution  required 
to  discern  them,  R,  is  equal  to, 

R  =  (area)"1/2.  (10) 

The  relation  of  equation  10  is  only  appropriate  where  the  two  patches 
are  completely  independent.  In  any  real  system  involving  lenses, 
chemistry,  anc  emulsions,  this  relation  is  only  valid  at  very  low 
spatial  frequencies. 

If  the  average  flux  per  area  is 

EXB  =  np/area, 

where  np  is  the  absolute  number  of  photons.  Equation  10  is  multiplied 


by  unity, 


Usi  ng 


R  =  (  1/  (EXB  )1/2)  (EXB)1/2  /  (area)1/2.  (12) 


Mmin/k  =  (EXE)'1/2,  (13) 


Equation  10  can  be  rewritten  as 


R  =  Mmin  (EXB)l/<:  /  k.  (14) 
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The  implicit  limitation  of  these  three  factors,  minimum  discernable 
modulation,  average  exposure,  and  resolution,  is  due  to  photon  noise 

C 

and  can  be  seen  in  figure  2  .  The  eye  normally  sees  a  time  averaged 

or  integrated  level  of  the  light.  In  low  light  photography,  the 
interval  of  integration  and/or  number  of  available  photons  is  reduced 
and  the  random  nature  of  the  photons  is  clearly  seen.  Photon  noise  is 
the  ultimate  lower  limit  of  the  image's  ability  to  convey  information 
and  experimentally  is  inseparable  from  grain  noise. 


(a)  3  x  103  photons  (b)  1.2  x  10^  photons  (c)  S.3  x  lO^photons 
(d)  7.6  x  105  photons  (e)  3.6  x  10^  photons  (f)  2.8  x  107  photons 
Figure  2,  The  effect  of  increasing  photon  exposure  on  image  quality 

This  derivation  assumes  the  exposure  distribution  will  be  Poisson 
in  nature  regardless  of  the  size  of  the  areas  considered  and  that  each 
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patch,  or  photon  detector  is  entirely  independent.  Those  assumptions 
aside,  in  general  the  ability  to  resolve  an  object  presented  in  the  KTF 
of  a  film  system  should  improve  with  increasing  exposure.  The  relative 
improvement  of  a  CPA  exposed  emulsion  is  relative  to  the  CPA 
contribution  to  the  average  exposure. 

There  is  one  final  consideration  to  apply  in  postulating  the  effect 
of  CPA  on  KTF,  Within  a  linear  system,  a  Fourier  transform  of  a  sum  is 
equal  to  the  sum  of  the  transforms.  This  rule  is  appropriate  to  a 
consideration  of  CPA  because” the  resulting  image  is  the  sum  of  tne  lens 
induced  spread  and  the  generalized  flare,  or  irradiance  from  the  CPA 
device.  The  line  spread  function  of  a  system  avfected  by  flare  can  be 
shown  as: 


1  c 

fig  3,  Spread  function  affected  by  flare 
The  resulting  spread  function  is  the  sum  of  the  flare  and  the  lens 
induced  spread.  T)  is  equal  the  ratio  of  the  energy  in  the  flare  to  the 
energy  in  the  optic  spread.  The  result  of  increasing  values  of  *1)  on 
the  system  KTF  is  shown  below10: 


fig  <,  KTF  of  a  system  affected  by  flare10 
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IMAGE  QUALITY  l€ASUREMENTS 

It  is  now  appropriate  to  focus  attention  upon  some  previous  image 

quality  studies  concerning  CPA  and  discuss  some  measures  of  image 

quality  that  are  available  for  this  study. 

Historically,  studies  of  CPA  effect  have  focussed  upon  the 

additional  film  speed  obtained.  Two  image  quality  studies  have  been 

performed.  One  used  resolution  targets  with  the  result  that  CPA  didn't 

affect  resolution  in  regions  where  normal  exposure  also  rendered  an 
4 

image.  A  study  at  Rochester  Institute  of  Technology  by  B.P..  Desai 
studied  the  effect  of  CPA  on  detective  quantum  efficiency  (DQE) 

c 

relative  to  normal  exposures.-  DQE  will  be  derived  in  the  next 

section.  The  result  of  that  study  was  that  CPA  produced  a  better  image, 
as  measured  by  a  higher  DQE  in  the  exposure  regions  where  a  normally 
exposed  image  was  only  marginally  possible. 

Other  measures  exist  to  evaluate  image  quality  which  will  be 
defined  or  derived.  Because  of  the  historic  significance  of  DQE,  that 
measure  will  be  derived  and  used  in  the  study.  A  measure  of  modulation 
transfer  function  known  as  the  second  moment  approximation  will  be 
derived.  Weiner  spectrum  is  a  measure  of  the  emulsion  noise  as  a 

function  of  spatial  frequency.  It  will  be  defined  and  incorporated  into 
the  derivation  of  the  final  measure,  information  content.  These  four 
measures:  DQE,  KTF,  Weiner  spectrum,  and  information  content  will 

yield  a  more  comprehensive  evaluation  of  CPA  than  previously  performed. 

DETECTIVE  QUANTUM  EFFICIENCY 

One  measure  of  image  quality  that  incorporates  contrast  and  grain 
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noise  as  a  function  of  exposure  is  called  DQE.  An  equation  for  DQE  has 
been  derived  from  a  theoretical  consideration  of  photon  receptors.  This 
derivation  is  excerpted  from  Dainty  and  Shaw,  Image  Science^.  A 
number  of  assumptions  are  made  that  detract  from  absolute  rigor  in 
favor  of  reasonable  application  and  these  assumptions  will  be  addressed 
at  the  end  of  the  derivation. 

Given  an  array  of  equal  sized  (area=a)  independent  photo  receptors 
irradiated  with  an  average  of  q  photons  per  detector,  the  probability 
that  a  single  detector  will  receive  r  photons  is, 

P(r)  =  qr  e'q/ri.  (15) 

For  the  purposes  of  this  derivation  it  is  assumed  the  detector  is  100£ 
efficient  in  photon  absorption  between  a  threshold  level,  T,  and  a 
fixed  saturation  level,  S.  Below  T  and  above  S  the  detector  will  be  0 % 
efficient.  While  this  detector  would  be  impossible  in  reality,  it 
roughly  approximates  a  photographic  emulsion.  A  photographic  emulsion 
requires  a  minimum  level  of  irradiance  to  achieve  a  density  above  the 
base  +  fog  level.  The  base  +  fog  level  is  that  level  of  density 
inherent  in  the  film  that  would  be  seen  if  the  totally  unexposed  film 
were  developed.  Conversely,  there  is  a  maximum  film  density  that  can  be 
achieved,  after  which  the  addition  of  more  photons  does  not  yield  any 
increase  in  density. 

In  any  exposure  the  number  of  detectors  that  will  saturate  is, 

CC 

nsd  =  2  qr  e_q/r:-  |ie) 

r=S 


! 


I 
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The  average  count  level,  cl,  is, 

S-l 

cl  =  V(r-T+1)  qr  e-fl/rl  +  Y  DL  qr  e“<?/r! .(17> 


where  DL  is  the  difference  between  saturation  and  threshold  level  plus 
one  or  DL  =  S-T+l.  The  equation  can  be  simplified  to. 


T-l 


cl  =  DL  (1  -  f ■)  e-S), 
T  T+l 


where  f-j  =  1  /DL (  qr/rl  +  ^Tqr7r!  +  ^Tqr/rl 

r=0  r*0  r«0 

S-l 

+  ../>V/r:)(18) 

r=0 

It  is  necessary  evaluate  the  gradient  or  change  in  detector  count  with 
change  in  input  photons,  which  is, 

g  =  d(cl  )/d( q)  =  DL  e'q(f]  -  dff^/dfq)).  (19) 

The  summations  simplify  to, 


S-l 

f2  *  f]  -  Ciff-J  )/d(  q)  =  1  /DL  V  qr/r ! .  (2 

r-T-1 


:0) 


Therefore 


where, 


g  =  DL  f?  e*q,  (21) 


EL-1 

f2  =  l/DL  ^  qr/r.' . 

The  gradient,  g,  will  have  a  value  of  0  until  q  exceeds  T  and  a  value 
of  unity  until  g  approaches  S,  the  saturation,  afterwards,  the  gradient 


will  fall  to  0  again. 

There  will  be  a  spatial  distribution  in  the  image  resulting  from 
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the  statistical  distribution  of  photons.  Of  interest  is  d(c1  )  or 
mean  square  fluctuations  in  counts.  This  can  be  calculated  from  the 
second  moment  of  the  distribution.  The  second  moment  about  the  mean 
level  is  equal,  tc  the  second  moment  about  the  origin  minus  the  square 
of  the  first  moment. 

c(02)  =  m2  =  m2'  -  cl2  (22) 

The  first  moment  is  simply  the  mean  count  level.  The  second  moment 
about  the  origin,  m2‘  is 


DL-1 


(23) 


The  summations  can  be  manipulated  to  yield 

m2‘  =  DL2(l-f3  e'q), 

where 


_  T  T+l 

f3  ■  l/DL2('^qr/r.'  +  3^qr/rl  + 


T-l 

v-1 


L~ j 
r=0 


C-i 

r=0 


r=0 


S-l 

+  (2DL-l)^qr/r:  (24). 

r*0 


Combining  equations, 

d(cl )2  =  DL 2 ( ( T - f 3  e'q)  -  (1  -  f]  e'q)2),  (25) 
where  f^  is  eq  24  and  f^  is  eq  16. 

2  2 

A  relation  of  variance  in  to  variance  out  or  d(  cl)  to  d  (q )  must 
be  defined.  The  noise  in  the  image  is  referred  back  to  exposure, 
d( cl )2/( d( cl )/d( q) )2  =  d(cl)2/  q2. 
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It  is  assumed  the  distributed  photons  follow  a  Poisson  distribution, 
therefore, 

d(q)2  =  q  (26). 

and  at  exposures  below  saturation,  the  ratio  of  signal  variance  in  to 
signal  variance  out  is  defined  as  DQE,  where, 

DQE  =  d(q)2/(d(ri)2/g2).  (27) 

Assuming  a  Poisson  relation  between  d(q)  and  q,  DQE  now  becomes, 

DQE  =  q  g/d(cl )2. 

The  ratio  DQE  will  always  be  less  than  1.0  because  of  entropy 
considerations.  DQE  can  be  shown  to  be  the  ratio  of  signal/noise  into 
the  imaging  system  to  signal/noise  out.  If, 

SNR  =  q  /Jq  =  .  (28) 

Therefore,  SNR2p  =  q.  Ey  the  same  relation,  the  noise  in  the  image 

is, 

SNRout  =  <(cV[d(q)/dUl)]2  =  d(cl  )2/g2.  (29) 

Therefore, 

DQE  =  d(q)2  g2/d( cl )2  =  q  g2/d(cT )2,  (30) 

which  is  the  ratio  of  the  squares  of  SNRout  to  SNRin. 

The  factors  of  DQE:  gradient,  g,  exposure  in  photons,  q,  and 
output  variance  d(c l)c  will  be  defined  in  terns  of  the  photographic 
image. 

In  photographic  systems,  the  image  is  evaluated  in  terms  of 
density,  which  is  defined  as, 

D  =  log10(Io/It)  =  -log10  T,  (31) 
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where. 


Io  =  irradiance  out 
It  =  irradiance  in 
T  =  transmission  =  Io/It 

If  the  light  subtracted  from  It  is  proportional  to  the  mean  count 
level,  cl,  then, 

Io  =  It  -  cl  b  It.  (32) 

where  b  is  a  fixed  proportion  of  the  incident  irradiance.  By 
substitution. 


D  =  -log10(l  -  cl  b),  (33) 

and  using  the  mean  count  level, 

D  =  log10(l  -  b  DL( 1  -  f1  e'q)).  (34) 

In  a  photographic  system,  the  gradient  is  defined  in  terms  of  density 
or  d(D)/d( f(q) ). 


Gamma  =  d(  D  )/d(  1  og-j  Qq ) .  (35) 

Because, 

d(log10q)  =  log10e  d(q)/q,  (36) 

Gamma  =  (q/log10e)(d(D)/d(q) ),  (37) 

and 


d(D)/d(q)  =  Gamma  log^e/q.  (38) 

The  scanning  aperture,  A,  covers  N  of  small  area  a,  detectors,  A  = 
Na.  Density  variance  is  now  being  measured  as  a  function  of  this 

O 

scanning  aperture  size,  which  is  referred  to  as  d(D  )  .  Likewise 

d 

exposure  is  considered  as  a  function  of  the  same  scanning  area,  giving 
qa:  da  =  q/A.  The  system  output  noise  is  now, 

d(D)2[d(q)/d(D)]  =  d(D)2[qa/Gamma  log^e]2.  (39) 

DQE  now  becomes, 

DQE  =  C(log10e)2Gamma2]/[qa  d(Dfl)2].  (40) 
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To  evaluate  d(P  )  ,  it  is  known  that, 

O 

D  =  -log1Q(l  -  b  cl).  (41) 

If  the  change  in  density  with  respect  to  change  in  count  remains 
constant  with  area, 

d(D)/d(cl )=  b  log10e/(l-b  cl)  =  d( )/d( cl a ) .  (42) 

The  finite  fluctuations  must  be  small  enough  to  make  this 
differential  equation  predictable.  This  can  be  done  by  making  the 
scanning  area,  A,  large  enough  to  encompass  a  large  number  of 
individual  detectors.  The  mean  square  fluctuation  measure  of  A  area 
containing  N  detectors  will  be  related  to  the  mean  square  fluctuation 
when  measured  for  individual  detectors  by, 

d(cl  )2  «  1/N(d(cl2)).  (42) 

a 

Considering  than  A  =  Na  and  combining  equations, 

A  d(Dfl)2  =  [a  (b  log10e)2  d(cl)2]  /  [1-b  cl].  (44) 

By  substitution, 

DQE  =  (logine)2  Gamma2 

qa  A  d(Da)2  .  (45) 

This  is  a  very  effective  tool  for  evaluating  an  emulsion's 

—  2 

capacity.  A,  d(Da)  ,  Gamma,  and  qa  are  easily  measured.  As 

«  d 

previously  indicated,  DOE  measurement  has  been  employed  in  the  past  in 
a  comparison  of  normal  vs  CPA  exposure.  A  Poisson  distribution  viewed 
through  a  finite  aperture  is  not  rigorously  Poisson  and  except  in  a 
highly  idealized  case,  detectors  do  not  possess  the  independence 
assumed  in  this  derivation.  Because  of  these  assumptions,  the  material 
is  practical  but  not  totally  rigorous. 


MODULATION  TRANSFER  FUNCTION 

In  this  section  the  basis  for  MTF  will  be  defined  and  a  simple 
approximation  derived  that  will  be  used  in  the  experiment.  Modulation 
transfer  function  (MTF)  is  a  tool  for  evaluating  the  relative  ability 
of  a  system  to  transfer  signal  information  as  a  function  of  spatial 
frequency.  To  understand  MTF,  a  one  dimensional  model  of  the  imaging 
system  is  assumed,  a  line  function  which  approximates  a  deltafunction 
in  one  dimension  is  input.  The  output  of  the  imaging  system  is  the  line 
spread  function,  l(x),  as  shown  below: 


figure  i.  Line  spread  function  output  of  an  imaging  system 

The  optical  transfer  function  (OTF)  is  the  Fourier  transform  of  the 
line  spread  function  or, 

CTF(f)  =  f llxle'1  27Tf  xdx.  (46) 

•L  OC 

and  the  MTF  is  the  modulus  of  the  OTF  or, 

MTF(f)  =  JoTF ( f )J  .  (47) 

The  one  dimensional  model  represents  only  a  slight  simplification.  If  a 
two  dimensional  delta  function  is  applied  to  the  system,  the  resulting 
impulse  response  function  is  the  point  spread  function,  P(x,y).  In  the 
absence  of  asymmetric  optical  aberrations,  P(x,y)  can  be  assumed  to  be 
radially  symmetric  about  its  peek.  The  line  spread  function  can  be 


l(x)  =  d(e(x))/dx  .  (49) 


found  by  integrating  the  point  spread  function  over  all  y, 

Kx)  =  P(x,y)dy.  (48) 

The  line  spread  and  point  spread  function  are  difficult  to  generate 
in  the  laboratory.  For  ease,  the  edge  function,  e(x)  is  used,  created 
by  a  light/dark  transition.  This  is  related  to  the  line  spread  function 


The  method  of  calculating  MTF  to  be  used  in  this  study  is  the 

second  moment  approximation.  Its  benefits  are  simplicity  and  accuracy 
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in  predicting  MTF  tc  about  0.30  restive  modulation.  To  calculate 
the  second  moment  MTF  approximation,  e^T’^x  is  series  expanded  , 


-i2T7fx 


=  t  (-i  2  TTf  x)n/n! 
n=u 


Fitting  this  into  equation  46  gives: 

KTF(f)  =  J  l(x)  ] 


f  Kx)  f;(-27Tfx)n/n:  dx 
-oc  n*0 


term.s 


;l  CO  00 

l(x)  dx  -  i  2  TTf  /  Kx)  x  dx 

oc  -cr 

-  2  T TZ  fZJ*l(x)  x2  dx  - 


l(x)  x  dx  -  higher  order 


Convention  dictates  that  the  area  of  Kx)  be  normalized,  therefore 

f l(x)  dx  =  1.0.  (53) 

th  r*  N 

The  Frr  moment  of  a  function  is  given  by  J  Kx)  x  dx  (54) 

-CX 

By  proper  choice  of  a  coordinate  system  the  first  moment,  can  be 
set  to  zero  and  the  second  moment,  M^  properly  computed.  The  equation 


for  MTF  now  becomes, 


MTF(f)  =  1  -  2TT2  f2  (55) 


WEINER  SPECTRUM 


The  Weiner  spectrum  is  a  measure  of  the  emulsion  noise  as  a 
function  of  spatial  frequency.  It  is  a  specialized  application  of  the 
concept  of  power  spectrum  which  may  apply  to  either  the  signal  or  the 
noise  in  an  emulsion.  The  Weiner-Khintchin  theorem  relates  the  image 
variance  i n  the  spatial  domain  to  the  spatial  frequency  domain.  The 

r 

Weiner-Khintchin  theorem  states  that  C(x)  and  W(f)  are  Fourier 
transform  pairs  where  C(x)  is  the  auto-covariance  of  any  function,  f ( x ) 


C(x) 


■P 

00 


(f{a)-u)(f  (a  +  x)  -  u)da  .  (56) 


where  f  (a)  is  the  complex  conjugate  of  f ( a )  and  u  is  the  mean  value 
of  the  process  or  DC  level.  The  transform  pairs  are. 


M  CC 

C(x)  =J  wme1  2TTf  x  df,  (57) 


and 


f  °° 

W(f)  = J  C(x)e_i  2TTf  x  dx.  (58) 

OC 

In  practice  the  Weiner  spectrum  is  measured  by  sampling  the  grain  noise 
with  an  elongated  slit.  The  sampling  interval,  dx,  of  the  slit  in 
conjunction  with  the  slit  width,  determines  the  maximum  detectable 
frequency  .  The  highest  detectable  frequenty,  fmax,  is  determined  by 

fmax  =  (2dx)”^ . '  (59) 

The  slit  width  is  convolved  with  the  grain  noise  during  the  sampling 
process.  The  sampling  slit  width  is  set  so  that  the  sine  function  that 
results  from  the  transform  of  the  rectangle  function  will  experience 
its  first  zero  at  fmax.  This  is  done  by  setting  the  slit  width  equal  to 
2dx.  The  purpose  of  this  selected  slit  width  is  to  prevent  alaising 
during  the  Fourier  transform  process.  Alaising  is  the  process  whereby 
high  frequency  power  becomes  confused  with  low  frequency  power.  In  any 


convolution,  of  which  auto-covariance  is  a  specialized  form,  the  width 
of  the  resulting  function  on  the  axis  is  equal  to  the  sum  of  the  axis 
widths  of  the  functions  being  convolved.  To  allow  for  this  increase  in 
required  space,  trailing  zeros  are  added  to  the  samples  taken. 
Specifically,  if  a  128  point  transform  is  to  be  performed,  64  points 
will  be  sampled  data  points,  64  points  will  be  zeros. 

Two  more  processes  are  applied  to  the  data  array.  The  first  is  to 
determine  the  mean  level  and  subtract  that  value  from  all  the  data 
points.  Secondly,  the  Fourier  transform  requires  that  the  function  to 
be  transformed  be  finitely  integratable.  The  sampled  noise  function 
doesn't  satisfy  that  condition.  To  achieve  a  finitely  integratable 
function,  a  "window"  is  applied  to  the  function  to  drive  the  function 
to  zero  within  the  sampled  space.  The  “window"  chosen  for  this 
experiment  is  the  Bartlett,  which  superimposes  a  triangle  function  upon 
the  sampled  data.  In  the  end  the  sampled  data  set  will  look  like: 

► - 4f§ 

0  N/4  N/2  ■ 

figure  6,  Bartlett  window  applied  to  Noise  data 

The  transform  from  the  spatial  domain  to  frequency  is  done  using 
the  Tukey-Cooley  Fast  Fourier  Transform  algorithm.^  A  real  and 
imaginary  value  results  at  each  data  point.  The  power  is  determined  by 
taking  the  square  of  the  modulus  of  the  complex  value.  The  value  of  the 
spacing  between  the  data  points  in  frequency  is  determined  by  the 
relation, 


N  df  dx  =  1.0,  (60) 
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or, 

df=  (N  dx)"1. 

Where  N  is  the  number  of  points  in  the  transform,  dx  is  the  sampling 

interval  and  df  is  the  data  interval  in  frequency. 

Because  noise  is  a  random  process,  the  correct  application  of 

statistics  is  critical  to  insure  the  noise  power  spectrum  is  correctly 

18 

measured.  The  standard  error  in  measuring  Weiner  Spectrum  is  roughly 

S.E.m  =  (X  df)'1/2.  (61) 

Where  X  is  the  total  scan  distance  and  df  is  the  data  interval  in 
frequency.  The  values  chosen  will  be  discussed  in  the  experimental 

section.  The  concept  of  Weiner  spectrum  will  be  used  in  the  final 
measure  of  image  qualit>,  information  content. 

INFORMATION  CONTENT 

12 

Information  science  began  with  the  initial  work  of  C.  Shannon 

and  gives  considerable  insights  into  the  content  and  capacity  of  an 
information  system.  Primarily,  information  is  a  function  of  the 
unexpected-ness  of  an  event.  Secondly,  there  is  an  additivity 
requirement  of  information. 


figure  7,  Additivity  requirement  of  Information 
If  thp  pvpnts  arp  i  nrlr-npndpnt.  SO  that  the  Drobabilitv  Of  both  is 


lltflll] 
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and  the  information  is  additive, 

I (A,B)  =  1(A)  +  1(B).  (63) 

The  combination  of  requirements  gives, 

I(A,B)  =  f { P (A )  P(B) ).  (64) 

The  only  function  with  this  property  is  the  logarithm,  therefore, 

KA.B)  =  k  logb(P(A)  +  P(B) ).  (65) 

For  I  increasing  as  P(A,B)  decreases,  k=-l  and  viewing  all  information 
as  ultimately  binary,  t*2.  Finally  our  expression  becomes, 

KA.B)  =  -log2(P(A)  +  P(B) ).  (66) 

If  an  image  with  Na  photo  detector  elements/unit  area  is  considered, 
each  of  which  may  exist  at  M  number  of  different  levels.  Further,  all 
levels  are  equally  probable  and  adjacent  detectors  are  independent. 
Then  the  information  capacity  in  the  image  will  be, 

I  =  -log2(P(given  array  of  pixel  values)),  (67) 

I  =  -log2(K  ‘  Na), 

=  +  Na  log2(M),  (68) 

This  is  only  a  theoretical  limit  of  information  and  in  no  way  is  an 
image  quality  descriptor. 

An  image  quality  descriptor  that  considers  both  noise  and 
moculation  across  the  dimension  of  spatial  frequency  is  information 
coincident.  Ar.  expression  for  information  content  can  be  derived  by 
considering  a  stochastic  recording  media  where, 
np  =  #  of  pixel  values  possible 

• 

S  =  #  of  stochastic  particles  per  pixel,  ie  AgX  crystals 

pi  =  probability  of  observing  a  pixel  value  including  the 
true  signal  and  the  underlying  stochastic  process. 
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fi  =  the  probability  of  an  observed  value  given  the  true  value 


If, 


POP  =  probability  of  an  observed  picture. 

POP  =  number  of  ways  a  sequence  of  pixel  values  can  occur 

without  repetition  times  the  probability  of  the  observed 
sequence  occurring 

POP  =  degrees  of  freedom  in  the  process 

Then, 

POP  =  SI  f^1  f|P2....f^np 

(Spl):  (Sp2K . (Spnp):  .  (6S ) 

This  is  a  multinomial  distribution  which  can  be  better  understood  if 

the  signal  is  a  binary  process,  ie  np  =  2  (yes  or  no) 

POP  =  S.'  f^1  (1  -  fn ) (S  “  Sp1 1 

_(Spl)i  (S-SplK  (70) 

The  information  in  the  picture  is. 


I  *  lim  (S  -  «c  ) ( 1/S  lcg2(P0P) 

=  1  i it ( S  -  *  )(l/s  log2CS: )  -  ffl/S  log2(Spi) 

1*1 

np 

+  SP1  log2(fi)  (71) 

Using  Stirling's  approximation  for  fa.torials  of  large  X, 

X:  =  ^2  TT  x  (X/e)X  (72) 


gives, 


I  =  lim(S- )[1/S  log2 /2TTs"  +  log2S  -  log2e 

1/S  log2  72TTS  pi  -  ]T,  pi(1og2Spi  -  log2e) 

1*1 


np 


+  ^  pi  1 og2 ( f i 


)] 


(73) 
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Applying  the  limit  and  canceling  terms  yields, 

np  m 

I  =  ^pi  log2(pi)  +  ^Tpi  log2(fi)  .  (74) 

Information  content  can  related  to  capacity  by  assuming  a  perfect 

ecorcing  or  fi  *  1.0  for  all  i,  the  equation  reduces  to, 

n£ 

1  =  ^  pi  T c 9 2 ^ P 1"  ^  •  (7E) 

And  if  each  state  is  equally  probable, 

I  =  np  log2(2).  (76) 

which  is  the  same  as  previously  defined  with  the  number  of  levels,  M  =2 
Information  on  film  can  be  represented  by  the  summing  process  shown 
below: 

© — © - © 

(nos) 

figure  8,  Imaging  X  onto  Y  with  the  addition  of  Noise,  NOS 


Y  =  X  +  NOS  where  X  is  the  object  imaged,  and  NOS  is  the  total  system 
noi se. 

pi  =  P(y) 
fi  =  P ( Y/X ) 

which  is  the  conditional  probability  of  Y  given  X.  Information  becomes, 
rm  np 

I  =  ^pi  log2(pi)  +  ^pi  1  og2 ( f i )  (77) 

In  terms  of  expectation  of  signal  and  noise, 

np  np 

I  =  ^P(Y)  1  og2 ( 1  /P ( Y ) )  -^PfY)  log2(l/P(Y/X))  (78) 

I  =  I (Y )  -  KY/X)  .  (79) 


The  distribution  of  P(Y)  must  be  found  that  maximizes  the  integral. 


MAX  PROE  =  J  P(Y)  log2(l/P(Y))  dY.  (81) 

This  is  found  when  P(Y)  is  Gaussian  or, 

P(Y )  =  1  /  exp[  -(Y  -  u)1 2/  20-2],  (82) 

Applying  a  Gaussian  P(Y)  and  P ( Y/X >  to  equation  82  at  their  greatest 
values  gives 

I  =  1/2  log2(  cr  2/  cr  JQS)  (83) 

1/2  log2((  cr  2  +  cr  Jos)/  <r  j]os) 

.  (84) 

2 

The  variance,  cr  ,  in  the  signal  and  noise  can  be  related  to  its 
appropriate  power  spectrum  through  the  Weiner  Khintchin  theorem  which 
was  addressed  in  the  discussion  of  Weiner  spectrum.  It  can  be  shown 
that, 

cr2  =  C(0),  (85) 


further  from  the  central  limit  theorem: 

O'2  =  C(0)  =  f  wlf)  df. 


(86) 


oo 


Information  content  can  now  be  expressed  as. 


1  =  1/2  /* log2C(Ps(f)  +  Pn( f ) )  /  Pn(f)],  (87) 

oo 

where  Ps(f)  and  Pn(f)  represent  the  signal  and  noise  power 
respectively.  Because  a  photograph  is  a  two  dimensional  information 
form,  the  information  in  a  unit  area  is  represented  by, 

I  =  1/2  log2[  1  +  S( v,u)/M( v, u) 3  dv  du  .  (88) 


*  •> 


Assuming  the  image  is  isotropic  and  converting  to  polar  coordinates 


2  2  2 

from  cartesian  with  f  =  v  +  u  .  The  equation  for  information 
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content  becomes, 

rfe 

1  =TT  Jo  1092[  1  +  S(f)/K(f);!  f  df  »  (89) 

where  fc  is  the  highest  critical  frequency. 

1 3 

R.  Clark  Jones  has  written  about  information  content  in  films 
and  his  work  is  based  upon  the  fore-mentioned  theory  using  equation  89. 
He  determined  that, 

S(f)  =  P(f)  KTF(f).  (90) 

where  P(f)  is  the  power  spectrum  of  the  object  imaged.  In  his  study  he 
used  a  set  density  level  in  the  film.  The  KTF  was  calculated  from  the 
point  spread  function  of  the  emulsion. 

There  are  some  limitations  in  Clark's  work  that  must  be  considered 
in  applying  information  theory  to  a  photographic  system.  The  first  is 
the  limitation  of  Fourier  analysis  only  to  linear  systems.  While 
information  is  conveyed  in  linear  units  of  exposure,  this  information 
is  perceives  in  units  of  density  on  the  photographic  emulsion.  Density 
as  a  function  of  exposure  is  non-linear.  However,  photographic 
emulsions  can  be  proven  to  be  quasi-1 inear  over  short  intervals  or  a 
density  tc  energy  transformation  can  be  performed.  If  the  density  to 
energy  transformation  is  made,  and  computations  performed  on  the 
appropriate  energy  values,  the  linear  requirements  of  Fourier  analysis 
are  satisfied.  However,  using  energy  values  diminishes  the  value  of 
information  content  as  a  measure  because  density  is  the  visually 
perceived  unit.  Therefore,  density  units  should  be  used. 

The  second  problem  is  information  coding  on  tne  image.  Information 
content  reaches  a  peak  value  when  the  signal  is  optimally  coded  which 
occurs  when  each  grain  becomes  a  unit  of  information.  The  grain  can 
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then  be  represented  mathematically  as  a  delta  function,  whose  power 
spectrum  usee  in  equation  50  becomes  a  constant  of  unity  for  all 
frequency.  While  theoretically  desirable,  this  information  coding  would 
appear  as  random  noise.  Any  larger  image,  which  would  be  an  aggregation 
of  many  grains,  would  have  a  power  spectrum  less  than  unity  for  all 
frequency. 

For  this  study,  an  edge  will  be  used  as  the  object.  The  power  of  an 

_2 

edge  in  frequency,  Pe(f)  is  proportional  to  (f)  .  To  prevent  PeC f ) 

from  rising  asymmptotically  to  infinity  as  f  approaches  zero,  Pelf) 
will  be  modified  to  give. 

Pelf)  =  A  /  (B  +  f2).  (SI) 

2  14 

The  value  of  B  =  (C.Qlfc)  is  found  from  experimental  history 

The  auto-covariance  of  an  edge  does  not  exist  in  a  rigorous  sense.  To 

employ  Pelf)  it  must  be  subjected  to  some  scaling  relative  to  the  edge 

being  analyzed  to  fulfill  the  condition  of  the  central  limit  theorem  or, 

cr2  =  f Pelf)  df.  (92) 

•ice 

If  the  edge  is  considered  as  an  array  of  points  as  shown  below, 
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figure  S,  Edge  as  an  array  of  points 
Then  the  variance  cf  that  array  can  be  ^alculated  from, 


c ' 2  =  lim(X-oc)  1/2X  ^d(Dlx))2  dx  .  (93 

If  DD  is  the  total  edge  difference,  the  d(D(x))=  PD/2  and  d(D(x)) 


DD  /4.  Applying  the  limit  and  integrating,  the  variance  of  the  edge 


array  is  now  DD  /4.  Because  of  the  value  of  B  chosen,  if  A  is  set 


equal  to  DD*/4  then  Pelf)  will  be  properly  scaled. 
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SUMMARY 

The  mechanism  of  CPA  has  been  shown  as  a  photon  addition  process  to 
achieve  a  minimum  stable  latent  image  state  with  reduced  scene 
irradiance.  Some  image  quality  measures  have  been  defined:  contrast, 
noise,  and  modulation,  and  postulated  the  effect  of  CPA  upon  each  one. 
Some  laboratory  measures  relating  to  these  image  quality  measures  have 
been  defined  cr  derived,  specifically:  DQE,  MTF  by  the  second  moment 
approximation,  Weiner  spectrum,  and  information  content.  The 
experimental  section  to  fellow  will  outline  how  these  techniques  will 
be  specifically  applied  in  this  experiment. 
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EXPERI RENTAL 

The  purpose  of  this  experiment  is  to  evaluate  image  quality  as  a 
function  of  exposure  in  a  comparison  between  normally  exposed  Tri-X  Pan 
emulsion  anci  Tri-X  Pan  exposed  with  the  benefit  of  concurrent  photon 
amplification.  The  hypothesis  is  that  there  will  exist  a  distinct 
exposure  region  where  normal  exposure  will  not  yield  an  image  and  the 
CPA  treated  emulsion  will  yield  a  better  image  in  that  region.  As 
exposure  is  increased  to  a  region  where  normal  methods  will  yield  an 
image,  the  normal  exposure  will  improve  relative  to  the  CPA  exposed 
image  and  then  surpass  the  CPA  image  in  quality.  The  diminished 
contrast  anc  increased  noise  of  CPA  exposures  should  be  the  prime 
contributor  to  this  result. 

In  discussing  the  coneuct  of  the  experiment,  there  were  four 
distinctive  pursuits:  1. Building  the  apparatus,  2. choosing  the  film, 
determining  its  sensitonetry  and  exposing  the  samples,  3. performing  the 
microdensitometry,  4. performing  the  necessary  data  reduction  of  the 
scan  oata.  First,  a  discussion  of  building  the  apparatus. 


EXPERIMENTAL  APPARATUS 

To  perform  the  experiment,  a  camera  equipped  with  a  CPA  device, 
some  artwork  tc  image  and  a  mounting  system  to  contain  all  three 
principal  pieces  were  needed.  The  School  of  Photographies  Arts  ar.d 
Sciences  at  Rochester  Institute  of  Technology  provided  a  suitable  Speed 
Graphic  camera  tody  capable  of  4  x  5  inch  film  format  and  a  Fuji  lEOmm 
f/B.C  lens. 
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The  CPA  circuit  and  lights  were  designed  and  fabricated  as  shown  in 

figures  10  and  11  respectively.  The  positioning  of  the  light  emitting 

diodes  (LED)  was  critical  to  evenly  illuminate  the  film.  The  lights 
were  placed  at  the  corners  of  a  3  3/6  inch  square  and  elevated  above 
the  film  plane  so  that  the  irradiance  in  the  corner  was  equal  to  the 
irradiance  in  the  center  of  the  square.  The  equations  for  the 

irradiance  at  points  A  and  B  and  the  drawings  below  define  the 

necessary  light  equilibrium, 

A  =  4  cos3(atan(2.3865/X),  (94) 

B  =  1  +  2  cos3(atan(3.37£/X) )  +  cos3(atan(4. 773/X) )  (95) 


figure  12,  Positioning  of  Light  Emitting  Diodes 
The  equations  were  derived  with  the  assumption  that  the  light  emitting 
dicdes  would  act  as  a  point  source.  To  achieve  the  desired  radiometric 
balance,  A  was  set  equal  to  B  and  the  equation  solved  for  X.  The  result 
obtained  was  3.3416  inches.  The  lights  were  adjusted  to  that  height.  In 
film  tests,  this  was  not  perfect  because  the  LED's  didn't  behave  as  a 
point  source.  However,  after  some  adjustment,  there  was  a  central  area, 
about  2"  in  diameter,  (figure  11)  where  there  was  no  significant 
density  difference.  All  measurement  work  was  done  in  that  central  area. 


CPA  CIRCUI 


SPEED  GRAPHIC  CAMERA 

WITH  CPA  CIRCUIT 
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figure  13,  CPA  irraciiance  pattern 

The  lighting  was  accomplished  with  a  0SRAM  8022  bulb  operated  at  1C 
volts  3.7  amps.  Power  was  stabilized  with  a  Thorensen  AC2000  power 
supply  with  current  flow  monitored  with  a  Weston  ammeter.  The  light 
source  was  used  in  two  ways:  First,  to  trans-illuminate  a  radial  step 
tablet  used  to  establish  the  film  sensitometry.  Secondly,  it  was 
configured  to  give  reflected  light  to  the  edge  art  work  as  shown  in 
figure  E.  Ir radiance  was  measured  in  lux  at  the  film  plane  using  an 
E.G.&  G.  radiometer.  Placement  of  the  radiometer  head  was  standardized 
with  wooded  jig  that  fit  in  place  of  the  film  holder  in  the  camera 
back.  Using  the  integrating  function  of  the  radiometer,  the  shutter  was 
calibrated  at  1/15,  1/30,  1/CO,  1/125,  1/25C,  and  l/500th  of  a  second. 
Using  film  and  the  E.G.&  G.  as  radiometers,  this  lighting  system  was 
found  to  be  extremely  stable  and  predictable. 

The  next  requirement  was  something  to  image.  A  standard  Kodak 
radial  step  tablet  was  used  to  determine  the  sensitometry.  The  radial 
step  tablet  was  necessary  because  of  irradiance  fall -off  from  the 
center  of  the  film  plane  to  the  edges.  Density  was  measured  at  a 
consistent  distance  from  the  center  of  the  pie  shaped  density  pattern. 

An  edge  was  necessary  to  measure  MTF.  The  artwork  was  made  by 
laying  a  razor  blade  across  a  piece  of  Kodak  2415  film  and  exposing. 
The  film  was  then  placed  emulsion  side  facing  a  Kodak  90£  reflection 


apparatus  layoui 


SPEED  GRAPHIC  CAMERA 


card.  The  purpose  of  this  arrangement  was  to  give  even  illumination 
across  the  light  side  of  the  edge  to  allow  ample  room  to  make  the 
Weiner  spectrum  measurements. 

The  mounting  geometry  was  shown  in  figure  14.  The  optical  conjugate 
relation  was  set  at  a  magnification  of  1/2.  All  devices  were  mounted 
upon  available  optical  bench  hardware  attached  to  the  research  darkroom 
counter  with  tape.  An  alignment  jig  was  made  for  the  Speed  Graphic  to 
insure  optical  alignment  was  net  oisturbed  during  camera  operation. 

FILM,  SENSITCMETRY,  AND  EXPOSURE 

Tri-X  Par.  emulsion  in  a  120  format,  emulsion  number  5063  819,  was 
chosen  as  the  test  film  for  two  reasons:  First,  this  film  has  figured 
prominently  in  former  CPA  studies.  Second,  in  film,  MTF  is  a  function 
of  film  speed.  The  higher  the  film  speed,  the  poorer  the  film  FTF.  The 
MTF  of  the  film  was  considered  in  the  belief  that  the  poorer  the  MTF, 
the  more  likely  the  MTF  would  be  degraded  by  the  CPA  effect.  Tri-X  had 
the  highest  film  speed,  therefore  the  poorest  MTF,  of  the  readily 
available  emulsions.  A  number  of  manual  development  processes  were 
tried  and  rejected.  It  was  found  that  a  Versamat  machine  processor  gave 
very  good  consistency  within  a  batch,  ie  the  variation  in  density  at  a 
set  exposure  level  was  less  than  +  .02  density.  The  film  was  processed 
as  a  batch  with  leading  and  trailing  sensitometric  strips  at  76°,  2.E 
feet  per  minute  transport  using  Versamat  641  chemistry. 

The  exposure  level  of  the  CPA  circuit  had  to  be  set.  .The 
recommended  exposure  was  one  that  produced  a  base  +  fog  density  0.1C 
density  units  above  a  normally  exposed  emulsion.  Exposure  time  in  the 
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CPA  device  is  controlled  by  resistor  R1 ,  which  was  a  decade  resistance 
box.  Exposure  at  various  resistance  levels  and  density  achieved 
linearly  regressed  against  resistance.  From  this  analysis,  a  resistance 
level  was  chosen  with  the  final  resulting  CPA  base  +  fog  density  of 
0.11  above  the  normal  base  +  fog  level. 

The  final  sensitometric  curves  are  shown  in  the  results  section. 
Based  upon  that  information,  the  experiment  would  use  three  exposures 
exposed  with  both  normal  and  CPA  at  -2.07,  -1.82,  -1.57,  -1.46,  -1.24, 
and  -.70  log  exposure.  This  was  intended  to  emphasize  the  low  exposure 
region.  Three  additional  exposures  were  to  made  using  CPA  only  at  -2.30 
log  exposure  because  no  normally  exposed  image  was  expected.  In 
actually  taking  the  edge  images,  two  exposures  were  diverted  from  -0.70 
to  take  one  exposure  each  with  normal  and  CPA  at  -0.93  log  exposure. 

The  lens  was  KTF  tested  at  f/5.6  and  f/11  at  the  infinite/focal 
length  conjugates.  Based  upon  those  results,  the  best  KTF  of  the  lens 
appeared  to  be  half  way  between  the  two  test  conditions,  or  f/8.  In 
taking  the  edge  images  the  log  exposure  was  determined  and  an 
approximate  exposure  time  selected.  A  90£  reflection  card  similar  tc 
the  card  used  in  the  edge  artwork  was  placed  in  the  object  plane  to 
give  a  uniform  reflection.  An  irradiance  value  was  calculated  and 
corrected  for  transmission  losses  within  the  film  on  the  artwork.  Then 
the  light  was  moved  along  its  rail  until  the  desired  irradiance  was 
measured  in  the  film  plane.  The  radiometer  was  carefully  removed,  the 
film  back  installed,  and  three  exposures  each  with  and  without  CPA 
were  made  of  the  edge  artwork.  The  process  was  repeated  until  the 
series  was  complete.  Finally  the  irradiance  was  measured  in  the  image 
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plane  as  a  function  of  wavelength  using  the  E.G.4  G.  radiometer  with 

2 

its  monochrometer.  This  irradiance  was  converted  to  photons  per  cm 
and  integrated  from  38C  nm  to  650  nm.  The  wavelengths  of  integration 
were  chosen  based  on  the  lower  limit  of  the  lamp  and  the  upper  limit  of 
the  film.  Then  the  film  was  machine  processed  as  formerly 
described. Because  of  the  variables  involved,  the  whole  process  had  to 
be  completed  without  interruption.  Numerous  runs  were  made  incomplete 
due  to  various  difficulties.  In  the  end,  every  event  was  successfully 
completed  and  the  desirec  exposure  series  gathered. 


MICRGDENSITCNETRY 

The  microcensi tometry  was  performed  on  a  Perkin-Elmer  PCS  system 

tied  to  a  PDP-11  for  data  processing.  Normal  KTF  measurement  of  an  edge 
is  performed  with  the  assumption  of  a  "global"  edge,  meaning  an  edge 
that  is  effectively  infinite  in  length  with  resolution  becoming  a  one 
dimensional  problem.  Within  imagery,  a  "global"  edge  is  often 

inappropriate  as  the  object  to  be  resolved  will  be  two  dimensional  with 
comparable  sizes  in  the  orthogonal  directions.  To  evaluate  the  effect 
of  noise  on  measuring  these  non-global  edges,  a  10  x  25  micron  slit  was 
chosen.  Forty  exactly  adjacent  scans  of  the  edge  were  made  with  a 
sampling  interval  of  3  micron  for  1024  points.  Considering  the 
convolution  of  the  slit  with  the  data,  the  single  slit  gave  a  two 

dimensional  frequency  domain  of  40  x  100  cycles  per  millimeter  (cpm). 

The  positioning  accuracy  of  the  microdensitometer  used  was  such  that 
the  scans  could  be  averaged  to  give  a  variable  frequency  domain.  If  two 
scans  were  averaged,  the  frequency  domain  available  was  20  x  100  cpm 
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and  if  four  scans  were  averaged,  the  domain  was  10  x  100  cpm.  Finally 
if  all  40  scans  were  averaged,  the  domain  was  1  x  ICO  cpm.  The  averages 
mentioned  were  those  made  in  the  processing.  The  10,  20,  and  40  cpm 
frequency  domain  roughly  corresponds  to  the  extremes  and  center  of  the 
human  visual  response  measured  at  the  cornea.  The  average  of  the  forty 
scans  was  the  "global"  edge  and  was  used  to  compute  the  MTF  for  the 
information  content  computations. 

Weiner  spectrum  data  samples  were  taken  with  a  124  x  6  micron  slit 
at  a  sampling  interval  of  3  micron.  Data  was  taken  in  blocks  of  64  so 
that  zeros  could  be  added  to  fill  a  128  point  transform.  105  blocks  of 
data  were  taken  to  give  a  total  of  6670  data  points  per  exposure.  This 
gave  an  X  distance  of  data  equal  to  20. Cl  rmi  which  gave  a  standard 
error  of  13$. 

SOFTWARE 

Computer  software  was  written  in  Fortran  77  to  perform  the 
necessary  computations.  Three  principal  programs  were  operated  from  a 
command  file.  The  first  cf  these  performed  a  forty  scan  average  of  the 
edge  data.  300  points  at  each  end  of  the  scan  were  averaged  to 
determine  Dmax  and  Dmin.  The  scan  was  then  evaluated  to  determine  the 
exact  extent  of  the  edge  by  value  comparisons  between  data  points  and 
the  Dmax  and  Dmin  values.  The  array  number  for  the  beginning  and  end  of 
the  edge  as  well  and  the  edge  midpoint  and  Dmax  and  Dmin  were  passed 
via  a  data  file  to  the  second  program. 

The  second  program  took  the  Dmax,  Dmin,  and  edge  coordinates  and 
made  KTF  computations  using  the  second  moment  approximation  for 
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singular  scans  and  averages  of  two,  four,  and  forty  scans.  The  maximum 
number  of  ^  values  that  could  be  calculated  from  the  available  data 
were  40,  20,  10,  and  1  respectively.  The  edges  were  always  scanned  from 
light  to  dark.  If  a  negative  value  of  ^  was  computed,  this  was 
considered  a  bad  edge  and  that  value  was  rejected  from  the  statistical 
computations  of  mean  and  standard  deviation.  The  number  and  values  of 
good  ^2 ' s  and  their  summary  statistics  for  all  considerations  were 
output  to  printer.  The  fort>  scan  averaged  Mg  and  the  edge  density 
difference  were  output  to  screen  and  used  in  the  third  program. 

The  third  program  read  the  noise  data  in  blocks  of  64  and  computed 
the  density  variance  and  the  mean  value  of  the  64  data  points.  The  mean 
value  was  subtracted  from  each  data  point  and  the  array  doubled  in  size 
with  the  addition  of  64  trailing  zeros.  The  Bartlett  window  was  applied 
and  the  128  point  array  Fourier  transformed.  The  power  of  the  noise  was 
achieved  by  taking  the  square  of  the  modulus  of  the  transformed  output. 
This  process  was  repeated  ICE  times  with  ensemble  averaging  of  the 
resulting  Weiner  spectrum.  To  determine  if  the  frequency  content  of  the 
Weiner  spectrum  shifted  with  exposure,  the  first  moment  of  the  array 
was  computed  from  0  to  fmax.  Finally,  the  program  prompted  for 
contrast,  exposure  in  photons  per  cmS  averaged  Kg,  edge  density 
difference  and  sample  identification  information.  The  program  computed 
DQE  and  information  content  and  output  Weiner  Spectrum  and  all  values 
read  during  the  prompting  and  computed  by  the  program. 

All  software  was  evaluated  with  known  input  data  to  insure 
computational  accuracy.  The  programs  are  attached  to  the  thesis  as 
attachments  1,  2,  and  3. 
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RESULTS 

The  results  given  in  this  section  are  presented  in  the  order  they 
were  developed  in  the  experiment  an/or  derived  in  the  theoretical 

section.  For  ease,  all  the  figures  referenced  in  the  results  section 

are  grouped  together  at  the  end  of  the  thesis.  The  data  conveniently 
groups  into  three  major  blocks  leading  to  the  principal  measures  of 
this  study:  1 .Sensitometry,  contrast,  density  variance,  and  DQE 
2.MTF  3.  Edge  step  difference  and  information  content.  In  many  cases 
the  standard  deviation  of  the  data  was  an  order  of  magnitude  less  than 
the  mean.  In  those  cases,  only  the  mean  value  at  an  exposure  level  is 
presented.  When  appropriate,  regression  analysis  is  performed  and 
plotted  with  the  data. 

Many  of  the  results  followeG  the  hypothesis  very  closely.  When  it 
is  appropriate,  this  observation  will  be  made  with  the  data.  The 
exposure  point  where  the  preferential  method  changes  from  CPA  to  normal 
will  also  be  addressed.  A  detailed  discussion  of  the  results  that  do 

not  directly  follow  the  hypothesis  and  ar,  overall  summary  of  the  data 

and  the  conclusions  to  drawn  will  be  presented  in  the  following 
section.  First,  the  results  leading  to  DQE  will  be  presented. 

FILM  SENSITOMETRY  AND  DQE 

Figure  15  shows  the  sensitometry  for  the  film  process.  There  is 
approximately  1  stop  or  0.30  log  exposure  difference  between  the 
processes,  the  speed  points  are  annotated  on  the  figure.  The  film  speed 
difference  is  lower  than  previously  published  6  stop  advantage  for  CPA. 
The  difference  lies  in  the  present  use  of  a  quantitative  measure  of  the 
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speed  point  given  by  the  exposure  level  that  gives  0.10  density  above 
base  +  fog.  A  subjective  measure  was  used  to  evaluate  the  film  speed 
difference  in  the  referenced  article. 

The  sensitometric  curves  were  fit  by  least  squares  to  a  quartic 
equation.  The  first  derivative  of  the  two  equations  were  taken  with  the 
results  plotted  in  figure  16.  This  shows  the  available  contrast  at  each 
exposure  level  and  is  the  beginning  of  a  result  that  continues  through 
much  of  the  data.  CPA  possesses  a  distinctive  exposure  domain  where  it 
will  produce  a  better  image.  As  exposure  increases,  the  normally 
exposed  image  improves  and  then  surpasses  the  CPA  image.  For  contrast, 
the  cross-cver  point,  or  log  exposure  level  where  the  preferential 
exposure  changes  from  CPA  to  normal  is  -1.92 

Figure  17  is  a  plot  of  density  variance  vs  log  exposure.  This 
measure  of  noise  in  the  signal  shows  CPA  to  have  a  significant 
disadvantage  until  very  high  exposure.  Figure  18  shows  DQE  vs  log 
exposure.  The  DQE  relation  shown  in  equation  45  is  very  sensitive  to 
change  in  contrast  or  gamma.  The  advantage  normal  exposure  possesses  in 
gamma  and  d(D)  ,  where  normal  exposure  will  produce  an  image,  is 
immediately  apparent.  The  dashed  extension  of  the  normal  exposure  DQE 
line  is  extrapolated  data  based  upon  the  exposure  level  where  gamma 
would  reach  zero.  CPA  holds  an  advantage  over  normal  exposure  only  below 
-2.C7  log  exposure. 

MODULATION  TRANSFER  FUNCTION 

Modulation  transfer  function  was  the  most  intensively  studied 
aspect  of  this  experiment.  Two  distinct  data  sets  are  presented: 
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first,  MTF  calculated  with  a  global  edge.  Second,  MTF  calculated  with 
increasing  orthogonal  frequency  domain  and  greater  noise.  First  the  MTF 
data  calculated  with  the  global  edge. 

Figures  19  and  20  present  the  second  moment  calculated  for  normal 
and  CPA  exposures  respectively.  A  smaller  value  of  Mg  means  a  better 
MTF.  The  data  was  fit  to  a  cubic  equation  by  the  least  squares  method 
with  the  regression  equations  shown  below,  le=  log  exposure. 

Normal , 

K2  =  4.9299  +  10.22  le  +  9.90  le2  +  2.79  le3 
R2  =  .5749 


CPA, 

M2  =  -8.0702  -  22.96  le  -  14.94  le2  -2.96  le3 
R2  =  .3165 

Figure  21  shows  the  regression  lines,  without  the  data  points,  together 
on  the  same  chart.  A  complete  discussion  of  the  factors  impacting  the 
measurement  of  MTF  will  be  presented  in  the  discussion  section  to 
follow.  Recalling  that  lower  Mg  means  a  better  MTF,  CPA  does  not 
possess  a  clear  advantage  in  the  low  exposure  domain.  Normal  exposure 
is  clearly  preferential  in  the  high  exposure  domain  beginning  at  -1.70 
log  exposure. 

Figure  22  shows  Mg  v.  log  exposure  with  various  slit  lengths  for 
the  normal  process.  The  dimensions  of  the  orthogonal  frequency  domains 
were  discussed  in  the  experimental  section.  As  the  ratio  of  the 
orthogonal  frequency  domain  increases,  there  appears  to  be  less 
exposure  dependence  in  the  value  of  Mg.  Alternatively  stated,  Mg 
seems  to  improve  with  exposure  in  the  least  averaged  or  most  noisy 
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scans.  As  the  averaging  decreases  the  noise,  doesn't  seem  tc 
improve  with  exposure  as  markedly.  Figure  23  shows  the  standard 
deviation  in  the  measurements  in  figure  22.  Figure  24  shows  the 
ratio  of  "good  edges"  that  were  used  in  computing  the  data  for  figure 
22.  Figures  23  and  24,  shew  in  general,  the  confidence  to  perceive  an 
edge  improves  with  increasing  exposure  and  diminished  orthogonal 
frequency  domai n. 

Figures  25  to  27  show  the  equivalent  results  for  the  CPA  process  to 
figures  22  to  24  for  the  normal  process.  CPA  results  are  generally 
similar  to  these  of  normal  exposure  with  an  exception.  The  ratio  of 
good  edges  does  not  follow  the  same  pattern  as  normal  exposure.  For 
short  slit  lengths,  the  ratio  of  good  edges  reaches  a  peak  less  than 
1.00  and  falls  of  slightly  with  increasing  exposure. 


INFORMATION  CONTENT 

Figure  28  shows  the  edge  difference  between  the  macro  maximum  and 
minimum  density  across  the  edges.  This  value  is  used  to  compute  the 
power  spectrum  of  the  edge  as  shown  in  the  theory  introduction.  The 
domains  are  clear  where  CPA  and  normal  exposure  are  advantageous  with 
the  cross-over  occurring  at  -1.46  log  exposure.  Recall  the  subjective 
six  stops  difference  between  CPA  and  normal  exposure  vice  the 
quantitative  one  stop  difference  used  in  this  experiment.  The 
subjective  measure  may  be  more  appropriate  considering  the  shallow 
slope  of  the  CPA  contrast  and  edge  difference  curves.  The  CPA  process 
may  give  a  acceptable  image  well  below  a  quantitatively  determined 
speed  point  making  a  six  stop  difference  very  reasonable.  This  effect 
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was  not  considered  in  establishing  the  sample  collection  scheme. 

Figures  25  and  30  show  the  information  content  vs  log  exposure  for 
normal  and  CPA  exposure  respectively.  The  information  content  was 
linearly  regressed  against  log  exposure  with  the  equation  values 
presented  below, 

Normal 

I.C.  =  387.51  +  151.98  le 
R2  =  .7945 
CPA 

I.C.  =  236.28  +  65.78  le 
R2  =  .6219 

In  the  comparison  of  these  two  charts  there  are  two  distinctive  domains 
wherein  each  technique  has  an  advantage.  The  cross-over  point  for  this 
measure  is  -1.84  log  exposure. 

Conspicuous  by  its  absence  is  data  concerning  the  Weiner  spectrum. 
The  first  moment  of  the  Weiner  spectrum  was  computed  and  found  not  to 
change  as  a  function  of  exposure.  From  the  central  limit  theorem  the 
integration  of  the  Weiner  spectrum  over  all  frequency  is  equal  to  the 
variance.  Further,  from  figure  17  the  density  variance  increases 
monotonically  with  exposure  at  a  rate  dependent  on  whether  normal  or 
CPA  exposure  is  used.  For  completeness,  figure  31  shows  a  typical 
Weiner  spectrum.  The  impact  of  the  Weiner  spectrum  on  the  information 
content  equation  is  ore  of  an  increasing  denominator  factor  with 
increasing  exposure.  The  power  spectrum  of  the  edge  is  generally 
increasing  as  a  numerator  factor  with  increased  exposure.  The  last 
numerator  factor,  MTF  in  some  domains  of  exposure  is  also  generally 
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increasing  with  exposure.  The  result  is  a  competitive  rate  problem. 
Judging  from  the  information  content  results  the  numerator  factors, 
edge  power  spectrum  and  MTF,  increase  faster  than  the  noise  as  a 
function  of  exposure. 

A  mere  complete  discussion  of  these  results  and  conclusions  to  be 
drawn  from  them  will  follow  in  the  next  section. 
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DISCUSSION 


INTRODUCTION 

This  discussion  section  will  begin  with  a  re-statement  of  the 
hypothesis.  There  is  a  distinct  region  in  which  CPA  will  provide  an 
image,  and  normal  exposure  will  not.  Within  that  region,  CPA  by  default 
provides  the  better  image.  There  exists  a  region  in  which  both  CPA  and 
normal  exposure  technique  will  provide  an  image.  As  the  exposure  level 
increases  ir.  this  region,  the  quality  of  the  normally  exposed  image 
v.'ill  improve  to  equal  and  then  finally  surpass  the  quality  of  the  CPA 
generated  image.  There  are  some  results  that  directly  support  that 
hypothesis.  Specifically  these  are  the  comparative  values  of  contrast, 
DQE,  edge  step  difference,  and  information  content  as  a  function  of 
exposure.  There  are  two  data  sets  within  the  results  section  that  are 
not  immediately  supportive  of  this  hypothesis.  These  are  density 
variance  and  bTTF  as  a  function  of  exposure.  For  brevity,  this  section 
will  only  address  the  data  non-comp! iant  with  the  hypothesis.  Secondly, 
the  nature  of  the  cross-over  point  will  be  discussed.  This  is  the 
exposure  level  at  which  the  preferential  exposure  means  changes  from 
CPA  to  normal.  First,  a  discussion  of  density  variance. 

DENSITY  VARIANCE 

As  previously  stated  density  variance  is  a  measure  of  noise.  This 
result  is  not  in  strict  compliance  with  the  hypothesis  but  is  a 
reasonable  result  because  CPA  has  added  noise  which  will  never  diminish 
below  the  base  level  achieved  on  the  film.  Therefore,  in  reviewing  this 
result,  shown  in  figure  17,  the  noise  level,  or  variance  increases  with 
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increasing  exposure  at  a  rate  dependent  upon  the  exposure  technique 
used.  The  results  also  show  the  noise  level  in  the  CPA  emulsion  is 
always  higher  than  the  normally  exposed  image,  but  the  difference 
decreases  with  increasing  exposure.  This  is  a  relatively  minor 
deviation  from  the  hypothesis  stated  and  does  give  interesting 
information  in  regard  to  the  additive  noise  given  by  the  CPA  effect. 

ktf 

KTF  does  not  fellow  the  hypothesis  directly  as  stated.  There  are  a 
number  of  reasons  for  this  and  include  both  errors  in  measurement  and 
real  changes  in  the  MTF.  Artificial  changes  in  MTF  result  from 
measurement  errors,  primarily  from  noise.  Noise  effect  is  inseparable 
from  the  measurement  but  one  can  postulate  the  effect  based  upon  the 
condition  under  which  the  KTF  measurement  was  made.  There  are  two 
possible  sources  of  real  change  of  the  MTF  which  will  be  discussed,  the 
first  of  these  is  the  change  in  the  perceived  point  spread  function  as 
a  function  of  exposure,  the  second  is  related  to  adjacency  effect. 

Noise  has  a  bi directional  effect  dependent  upon  the  power  of  noise 
relative  to  the  power  of  the  signal  and  the  variance  of  the  noise 
across  the  edge.  This  effect  can  either  help  or  hinder  the  KTF 
measurement.  Alternatively,  noise  will  either  give  an  inflated  measure 
or  deflated  measure  of  MTF.  In  the  case  of  an  inflated  measure  the 
measure  will  be  considered  inaccurate  and  not  representative  of  the 
true  KTF  but  rather  an  artifact  of  measurement.  In  the  case  where  noise 
would  hinder  the  MTF  measure,  this  a  realistic  measure  of  the 
operative  KTF  because  that  noise  effect  would  be  a  factor  in  any 
viewing  of  the  image. 
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Noise  can  be  a  help  or  enhance  the  measure  of  KTF  when  two  specific 
conditions  are  met.  One,  the  edge  being  scanned  is  relatively  low 
contrast  and  the  noise  is  statistically  stationary  across  that  edge. 
Second,  when  the  noise  power  is  significantly  less  than  the  power  of 
the  edge.  Under  these  conditions  the  noise  will  add  a  random  phase 
factor  with  the  average  result  being  an  improvement  in  the  measured 
KTF15. 

The  effective  noise  can  hinder  or  can  reduce  the  measured  MTF,  when 
two  conditions  mentioned  above  are  not  present.  Either  the  power  of  the 
noise  is  greater  than  the  power  of  the  edge  or  the  noise  is  not 
statistically  stationary  across  the  edge.  If  either  of  these  conditions 
are  present,  noise  will  hinder  high  frequency  signal,  thereby  reducing 
KTF. 

As  stated  in  the  introductory  paragraph  of  this  discussion  of  MTF, 
the  effective  noise  is  either  a  help  or  hindrance  and  becomes  a 
measurement  artifact  in  evaluating  MTF.  There  are  two  distinct  factors 
which  change  as  a  function  of  exposure  and  do  result  in  a  distinctive 
change  in  the  KTF.  The  first  of  these  is  the  change  in  the  point  spread 
function  as  a  function  of  exposure.  In  general,  the  line  spread 
function  of  an  entire  imaging  system  on  an  emulsion  will  take  the  form 
shown  below: 


figure  32,  Typical  Line  Spread  Function 
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This  form  will  vary  with  optical  aberrations  involved  and  chemical 
effects  induced  in  the  film.  Assuming  radial  symmetry,  the  key  point  is 
the  shape  of  the  point  spread  function  narrows  as  it  approaches  the 
peak  and  has  a  distinctive  slope  function  as  it  approaches  its  null. 
Within  a  photographic  emulsion  there  is  a  measurable  domain  given  by 
the  density  above  base  +  fog.  The  shape  of  the  point  spread  function  in 
a  rigorous  sense  will  not  change,  however,  the  measurable  domain  of 
that  point  spread  function  will  change  as  as  function  of  exposure. 

An  analogy  to  a  volcano  creating  an  island  is  very  appropriate.  The 
island  does  not  change  shape  in  an  absolute  sense.  However,  as  the 
volcano  generates  more  island  form,  so  that  the  island  emerges  above 
the  surface  of  the  the  water,  the  shape  of  the  volcano  becomes  broader 
to  those  viewing  from  above  the  water.  In  an  analogy  to  a  photographic 
system,  there  is  a  base  +  fog  level  on  the  film  which  equates  to  the 
water  in  the  example.  As  increased  exposure  elevates  the  point  spread 
function  of  the  imaging  system  above  this  base  +  fog  level,  the  visible 
domain  will  tend  to  broaden  until  the  first  null  is  reached.  After  the 
first  null  is  above  base  +  fog,  there  will  not  be  any  relative  increase 
in  the  width.  The  volcano  analogy  loses  some  of  its  value  because  the 
width  of  the  volcano  is  not  subject  to  any  finite  limitations.  The 
imaging  system  point  spread  function  does  have  a  limiting  size  based 
upon  diffraction  and  the  aberrations  of  the  system.  However  one  may 
generalize  that  KTF  will  diminish  slightly  with  increasing  exposure  due 
to  the  visible  or  measurable  domain  of  the  point  spread  function. 

There  is  an  ancillary  effect  within  photographic  emulsions  called 
20 

adjacency  effect.  The  adjacency  effect  is  created  by  chemical 


anisotropy  which  results  from  the  unique  nature  of  hydroquinone  based 
developing  solutions.  Where  the  edge  occurs  there  will  be  a  relative 
area  of  high  density  requiring  exhaustive  chemical  efforts  and  low 
density  requiring  relatively  little  chemical  effort.  The  nature  of 
hydroquinone  based  developers  is  such  that  where  there  is  increased 
developer  activity,  the  exhausted  developer  will  locally  change  the  pH 
of  the  development  solution  resulting  in  increased  developer  activity. 
The  apparent  effect  upon  the  emulsion  is  that  the  dark  get  darker, 
lights  get  lighter  at  the  boundary.  Agitation  will  reduce  this  effect. 
Machine  agitation  in  the  Versamat  should  normally  be  adequate,  however 
more  vigorous  methods  would  have  reduced  the  adjacency  effect  more.  The 
effect  is  shown  in  the  diagram  below: 

A  _  ftnax 

/ 
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figure  33,  Example  of  adjacency  effect 
The  net  result  is  that  this  will  tend  to  improve  MTF  with  increased 
adjacency  effect  with  increases  with  increasing  exposure. 

An  attempt  to  correlated  these  effects  to  the  results  is  very 
complex,  but  some  postulation  can  be  made.  In  viewing  results  figure 
IS,  at  low  exposure  the  noise  has  relatively  low  power  and  is 
statistically  stationary  across  the  edge.  The  result  is  an  improvement 
in  the  measure  of  KTF  due  to  this  random  phase  factor  given  by  the 
noise  under  these  conditions.  Secondly,  at  this  exposure  there  is  a 
small  apparent  point  spread  function,  again  resulting  in  improved  fTTF . 


As  the  exposure  increases,  three  factors  are  occurring,  one  the  noise 
loses  its  stationary  quality.  The  noise  power  also  becomes  greater.  The 
point  spread  measurable  domain  expands  and  finally,  the  adjacency 
effect  given  by  the  developer  improves.  Reviewing  the  results  we  can 
see  a  marked  rise  in  the  M^,  followed  by  a  rather  slow  decline  in 
K.2  or  improvement  in  MTF.  In  the  region  where  MTF  is  gradually 
improving,  the  predominant  effect  seems  to  be  the  adjacency  effect. 

In  reference  to  results  figure  22,  one  can  see  the  effect  of  noise, 
adjacency  effect,  and  point  spread  function  visible  domain  in 
measurements  where  the  power  of  the  noise  is  significantly  greater.  The 
inclusion  of  this  noise  destroys  the  available  high  frequency  signal 
resulting  in  a  diminished  MTF  measured.  As  the  noise  is  reduced  through 
greater  averaging  and  reouced  orthogonal  frequency  domain,  the  measured 
MTF  values  approach  those  measured  with  a  global  edge.  Therefore,  in 
the  presence  of  noise  the  ability  to  detect  a  two  dimensional  otject 
with  comparable  orthogonal  measurements  significantly  improves  with 
increasing  exposure.  As  the  ratio  of  the  object  dimensions  increase, 
the  exposure  dependence  decreases,  this  result  is  also  shown  in  results 
figure  24,  the  ratio  of  good  edges.  The  ratio  improves  with  increased 
averaging  and  increased  exposure.  Further,  the  relative  slopes  of  the 
lines  diminish  with  increased  averaging. 

Many  similar  arguments  apply  to  the  CPA  exposed  images  both  with 
the  global  edge  and  the  increased  noise  measurements.  Because  there  is 
a  unique  domain  in  which  only  CPA  will  give  an  image,  any  MTF  measured 
will  be  better  that  the  normally  exposed  image  measured  MTF.  The  M2 
increases  to  a  peak  level  and  then  begins  a  rapid  decline  to  converge 


upon  the  normally  exposed  value.  The  reason  for  the  shift  in  the  peak 
location  towards  higher  vs  lower  exposure  is  not  clear.  In  reference  to 
measurements  made  with  increased  noise  and  orthogonal  frequency  space, 
many  of  the  trends  seen  in  the  normally  exposed  images  are  apparent. 
There  is  an  exception  noted  in  the  results  section  concerning  the  ratio 
of  good  edges  not  converging  to  one.  System  noise  in  the  measured 
values  coupled  with  a  reduced  adjacency  effect  in  the  CPA  images  are 
possible  explanations. 

Considerable  emphasis  has  been  given  to  minute  differences  in  the 

measured  values  and  attempts  made  to  derive  conclusions  on  them  in 

relation  to  exposure.  Much  of  this  may  be  unfair  in  light  of  a 

statistical  analysis  cf  the  data.  In  both  normal  and  CPA  cubic 
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regressions  of  the  global  edge  data,  the  R  is  relatively  low, 
especially  for  CPA.  Linear  regressions  for  both  data  sets  included  a 
slope  of  zero  within  a  9C*  confidence  limit,  indicating  that  there  is 
no  exposure  dependence.  The  second  moment  approximation  for  KTF  is 
sound,  however  its  ability  to  accurately  depict  complex  changes  in  KTF 
given  with  the  addition  of  flare  may  be  questionable.  Alternatively,  a 
more  comprehensive  measure  of  MTF  would  not  as  easily  lent  itself  to  an 
analysis  relative  to  exposure.  With  reference  to  figure  2,  there  should 
be  some  measurable  improvement  of  MTF  with  exposure.  A  possible  future 
work  could  be  some  power  spectrum  analysis  of  edge  data  as  a  funtion  of 
exposure  to  more  clearly  determine  this  exposure  dependence  of  KTF.  To 
conclude  the  remarks  on  KTF,  there  may  be  an  exposure  dependence  of  MTF 
that  is  affected  by  the  variety  of  factors  mentioned.  Unfortunately, 
the  data  does  not  support  that  conclusion  with  the  greatest  strength. 


53 


CROSS-CVER  POINT 

There  was  not  a  clear  agreement  within  the  data  on  a  boundary 
exposure  where  the  preferential  exposure  method  changed.  The 
cross-overpoint  was  measurement  specific  and  there  was  a  range  in  which 
the  measurements  experienced  their  cross-over.  For  review,  they  are 
listed  below: 

Measurement  Cross-Over  Point  (log  exposure) 

Contrast  -1.92 

DQE  -2.07 

hi  2  -1.70 

Edge  Difference  -1.42 

Information  content  -1.84 

The  mean  value  of  the  cross  over  point  was  -1.79  log  exposure  with  a 
standard  deviation  of  0.246.  This  average  point  corresponds 
approximately  to  the  half  speed  point  or  0.15  log  exposure  below  the 
normal  exposure  speed  point.  Because  of  the  high  standard  deviation, 
one  would  infer  a  boundary  region  rather  than  a  point.  Further,  the  use 
of  CPA  in  an  exposure  region  below  the  highest  boundary  point  should 
yielo  the  preferred  image.  With  decreasing  exposure,  the  benefit  of  CPA 
relative  to  normal  exposure  will  be  enhanced. 
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CONCLUSIONS 

In  summary  three  distinct  conclusions  can  be  drawn  from  this  study 
which  represent  the  assertion  of  the  initial  hypothesis.  There  are: 

1.  There  is  a  low  exposure  domain  in  which  only  CPA  can  give  an 
image.  The  size  of  that  domain  is  a  function  of  the  means  of  measure. 
For  a  quantitative  measure,  the  size  is  modest,  only  one  stop  in  this 
study.  However  if  a  subjective  measure  is  employed,  that  CPA  domain  of 
preference  could  be  larger,  as  much  as  six  stops. 

2.  There  is  a  domain  where  both  CPA  and  normal  exposure  will 
produce  an  image  and  in  this  domain,  CPA  will  give  a  better  image  than 
normal  exposure.  The  benefit  of  CPA  diminishes  towards  the  high  end  of 
this  transition  domain. 

3.  Above  the  transition  domain,  normal  exposure  will  yield  a  better 
image  than  the  CPA  treated  image. 
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Attachment  1  code 


c 

C  CMTF.  FTN  PROGRAM  WRITTEN  BY  KONRAD  KERN 

C 

C  THIS  PROGRAM  MAKES  AN  INITIAL  LOOK  AT  THE  DATA  AND  DETERMINES  THE 

C  UPPER  AND  LOWER  BOUNDS  OF  THE  EDGE  SCAN. 

C 

C  PROGRAM  INPUTS 

C  40  EDGE  SCANS  OF  1024  POINTS  EACH 

C 

C  PROGRAM  OUTPUTS 

C  1.  DMAX  AND  DMIN  FOR  THE  40  SCAN  AVERAGED  EDGE 

C  2,  BEGINNING,  MID-POINT,  AND  END  COORDINATES  OF  THE  EDGE 

C 

C*****************#***«**********************#****##****#*#*##*****###*#*< 

R£AL*8  MRR ( 4  > , MR ( 40  > ,  EDGE ( 1 024  > .  AVE,  AVL,  AVH 
REAL*8  M2A, M22A, D£N< 10), LE< 10 > 

REAL*B  AEDCE< 1024),  Ml.  M2 
INTEGER  SAMPLE 
EXTERNAL  F0UR1 


CALL  ASSIGN<5,  'TI :  ') 

CALL  ASSIGN<2, 'RES  DAT ' ) 

CALL  ASS1GNO,  'EDGE  DAT') 

13  FORMAT (F5  2, X, F4. 2) 

11  FORMAT (2014) 

15  FORMAT (13) 

16  FORMAT (F5  2) 

12  FORMAT ( G18  10) 

C  BEGIN  READING  THE  SCANS  AND  COMPUTING  THE  SECOND  MOMENTS 
815  FORMAT (/, 5X,  'LOG  EXPOSURE* ',  F8  3) 

802  FORMAT (  /,  5X«  'SAMPLE  NUMBER*  '.  13,  10X,  'CONTRAST* ',  F5.  3 ) 

803  FORMAT(/, 5X,  'EXPOSURE*  ', G1 1  5, 5X,  'EXPOSURE  TYPE*  ',13) 

904  FORMAT  (  / ,  4X ,  '«'.3X,  'SECOND  MOMENT',  10X,  'ONE  SCAN') 

C  PERFORM  AVERAGE  OF  FORTY  SCANS  AND  COMPUTE  THE  M2 
REWIND  1 

WR ITE( 5. *) 'AVERAGE  OF  FORTY  SCANS' 

LS*0 
M2A*0  O 
M22A*0  O 

C  READ  ALL  FORTY  EDGE  SCANS  AND  AVERAGE 
DO  310  1*1, 1024 
310  AEDGE ( I ) *0.  0 

DO  321  1-1. 40 

CALL  MYF(EDGE,  1024,  1 ) 

DO  320  11*1, 1024 

320  AEDGE (II) *AEDGE (II) +EDGE (II) 

321  CONTINUE 

DO  331  11*1. 1024 

331  EDGE (II) *AEDGE (II) /40,  0 

C  DETERMINE  DMIN  OF  THE  EDGE 
AVE*0  0 

DO  325  11*1,300 
AVE=AVE*EDGE( 1 1 ) 

AVE-AVE/300.  0 


325 
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Attachment  1  code,  continued 


C  DETERMINE  DMAX  OF  THE  EDGE 
AVH-O.  0 

DO  326  I 1*723/ 1024 

326  AVH=AVH-*-EDGE ( 1 1 > 

AVH=AVH/300.  0 
AVL=AVE 

C  DETERMINE  THE  EDGE  STEP  DIFFERENCE 
DD=AVH-AVE 

WRITE (2. 12)  AVL, AVH,  DD 
AVE* < AVE+AVH > /2  0 

C  DETERMINE  THE  MID  POINT  OF  THE  EDGE 
DO  327  11*1, 1024 
K*II 

I F ( AVE-EDGE (II))  328,327,327 

327  CONTINUE 

328  CONTINUE 

C  DETERMINE  THE  END  COORDINATE  OF  THE  EDGE 
DO  350  11*1, 1024 

I F ( AVH-EDCE (1I)>355.355. 350 
350  CONTINUE 

355  LH-II+3 

C  DETERMINE  THE  BEGINNING  COORDINATE  OF  THE  EDGE 
DO  360  11*1,  LH 
LL*LH— I I 

1 F  ( EDGE  ( LL )  -AVL  )  365,365,360 
360  CONTINUE 

365  LL*LL-3 

DO  390  I*LL.LH 
J*I— LL+1 

390  EDGE ( J ) “EDGE ( I ) 

K=K-LL 

C  COMPUTE  THE  SECOND  MOMENT  OF  THE  EDGE 

CALL  SECMO < EDGE,  ftl.  M2,  AVL.  AVH,  J,  K> 

LS*1 

UK*LH-LL 

WRITE (2,  15)LL.  K,  LH 
WRITE (2,  12)Mi 

850  FORMAT <  1  X,  'LL*'.  14,  3X,  'SPREAD*'.  14,  3X.  'LH*'.  14) 

WRITE (5,  800)  LS.M1.M2 

800  FORMAT  (IX,  13,  3X,  'Ml*',  Cl  5  5,  3X,  'M2*',G15.  5) 

500  F0RMAT(3X, 015  5, 3X. 14) 

810  FORMAT ( / , 5X,  'EDGE  DENSITY  CHANGE*  '. CIO.  5, 5X,  'NOISE  SDEV*'.C10  5) 
999  STOP 

END 

C  SECOND  MOMENT  SUBROUTINE 

SUBROUTINE  SECMO (EDGE, Ml,  M2,  AVL.  AVH,  J,  K) 

R£AL»8  EDGE( 1024), ALINE (100).  MI.  M2 
REAL*8  AMI, AM2,  AVL,  AVH 
AM2*0.  0 
AM1*0.  O 

I F ( AVH- AVL  >  800 , 800 . 5 
5  DO  10  1*1, J 

1 0  EDGE ( I )*< EDCE ( I ) -AVL ) / <  AVH-AVL ) 
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Attachment  1  code,  continued 


C  CALCULATE  Ml 

DO  SO  1*2, J 
S=I-1 

SO  AMI  =AM1  ( EDGE  ( I )  -EDSE  (I— 1))*(  S*0.  003+0.  OOIS) 

WRITE ( 5, * )  'CALCULATED  Ml AMI 

C  CALCULATE  M2  THE  FIRST  TIME 
DO  60  1*2,  J 
S*I-1 

60  AM2*AM2+  < EDGE ( I ) —EDGE (I-1))*((0.  003*S-AMl+0.  001S>**2> 

GOTO  810 

800  M2*- 100 

GOTO  820 

810  M1*AM1 

M2*AM2 

820  RETURN 

END 

C  SUBROUTINE  TO  READ  THE  EDGE  SCANS 
SUBROUTINE  MYF(VAL,  N,  LF) 

INTEGER  EDGE (1024) 

REAL*B  VAL  < 1024 ) 

CHAR AC TER *80  HEAD ( 3 ) 

READ(LF. 900)  (HEAD ( I ) . 1*1 . 3 > 

900  FORMAT (A80) 

READ(LF, 901)  ( EDGE ( J ) ,  J* 1 ,  N ) 

DO  10  J-l.N 

10  VAL( J)*EDGE( JJ/1000.  O 

901  FORMAT ( 20 14) 

902  FORMAT (IX,  2014 ) 

RETURN 

END 


Repeat  40  Times 
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Attachment  2  code 


C*****«************************#********#**#*****#*******«*****#******4 

c 

C  MTF.FTN  PROGRAM  WRITTEN  BY  KONRAD  KERN 

C 

C  PROGRAM  INPUTS 

C  1.  DMAX  AND  DMIN  VALUES  FOR  THE  TOTAL  AVERAGED  EDGE 

C  COMPUTED  IN  FORMER  PROGRAM  * 

C  2,  EDGE  BEGINNING.  MID-POINT.  AND  END  ARRAY  COORDINATES 

C  3.  40  ARRAYS  OF  EDG£  SCAN  DATA.  1024  SAMPLE  POINTS  EACH 

C 

C  PROGRAM  OUTPUTS: 

C  1,  M2  VALUES  COMPUTED  USING  SINGLE  SCAN  VALUES.  STATISTICS 

C  2,  M2  VALUES  COMPUTED  USING  AVERAGES  OF  2  SCANS.  STATISTICS 

C  3.  M2  VALUES  COMPUTED  USING  AVERAGES  OF  4  SCANS.  STATISTICS 

C  4.  M2  VALUE  COMPUTED  USING  AVERAGE  OF  40  SCANS 

C 

C*******************************************«**************************4 


REAL*8  MRR  <  4 ) . MR ( 40 ) . EDGE ( 1024 ) . AVE. AVL.  AVH 
REAL*8  M2 A, M22A.  LS.  DEN ( 10  > .  LE ( 10) 

REAL*8  AEDGE< 1024).  Ml.  M2 
INTEGER  SAMPLE 
EXTERNAL  FOUR  1 
CALL  ASSIGN ( 5. 'TI  ') 

CALL  ASSIGN (4,  'RES  DAT') 

CALL  ASSIGN (2. 'MRES  DAT') 

CALL  ASSIGN ( 1 . 'EDGE  DAT ' ) 

13  FORMAT<F5  2,  X.  F4  2) 

11  FORMAT (2014 ) 

15  FOR MAT ( 13) 

16  FORMAT (F5  2) 

C  SAMPLE  IDENTIFICATION  INFORMATION 
WR I TE ( 5.  * )  'SAMPLE  NUMBER  ?' 

READ (5. 15)  SAMPLE 

WRITE  (5,*)  'NORMAL-1.  CPA-2' 

READ  <5. 15)  IDS 

WRITE (5. *> 'LOG  EXPOSURE-?' 

READ <5. 12)  EXL 

C  READ  IN  DATA  FROM  FORMER  PROGRAM 
READ( 4.  12  > AVL 
WRITE (5. * ) AVL 
READ ( 4, 12 ) AVH 
WRITE (5. *  > AVH 
READ( 4.  12  >DD 
WRITE (5. *>DD 
READ (4, 15)LL 
WRITE (5. *)LL 
READ (4, IS)  K 
WRITE (5. * )K 
READ  <  4, 1 5  >LH 
WRITE (5. * )LH 
READ  <  4.  12  >M1 
FORMAT (GIB  10) 
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Attachment  2  code,  continued 


C  BEGIN  READING  THE  SCANS  AND  COMPUTING  THE  SECOND  MOMENTS 
WRITE!2,  802)  SAMPLE.  CON 
WRITE (2. 815)  EXL 

813  FORMAT(/.  3X,  'LOG  EXPOSURE*  '•  F8.  3 ) 

802  FORMAT  (  /.  3X,  'SAMPLE  NUMBER-'.  13.  10X.  'CONTRAST* F3.  3) 
WRITE<2.  803)  QA.  IDS 

803  FORMAT(/. 3X,  'EXPOSURE-  Cl 1.  3. SX.  'EXPOSURE  TVPE*  '.13) 

WRITE (2. 904) 

904  FORMAT  ( /»  4X«  '#',3X.  'SECOND  MOMENT'.  10X.  'ONE  SCAN') 

C  SINGLE  SCAN  DATA  COMPUTATIONS  OF  M2 

24  LS*0 

M2A-0. O 

M22A—0, 0 

DO  20  1*1.40 

CALL  MYF(EDGE. 1024, 1) 

DO  37  II-LL. LH 
J-II-LL+1 

37  EDGE  <  J > -EDGE (II) 

CALL  SECMOCEDCE.  Ml,  M2.  AVL.  AVH.  J.  K> 

IF  (M2 )  20.20.29 

29  LS-LS+1 

M2A-M2A+M2 
M22A-M22A+M2*M2 
MR ( LS ) -M2 

WRITE(3,800)LS. Ml.  M2 

800  FORMAT! IX.  F4  1. 3X.  '  Ml-  '.  IX,  Cl  1.  3.  3X,  '  M2-  ',  IX,  Oil.  3) 

20  CONTINUE 

WRITE <2. 804)  (J.MR(J), J-l.LS) 

804  FORMAT! IX.  14. 3X. Cl  1  3) 

IFILS-2  0)  .<1.32.32 

32  D— SORT ! ABS ! LS*M22A-M2A*«2  > / <  LS* ! LS— 1 >  >  > 

M2A-M2A/LS 

31  WRITE! 2. 805)  M2 A. D 

805  FORMAT!/.  3X.  'AVERAGED  M2-  ' ,  01 1.  3,  3X.  'SDEV  -  ',©11.3) 

10  CONTINUE 

C  PERFORM  AVERAGE  OF  TWO  SCANS  AND  RECOMPUTE  THE  M2 
REWIND  1 
WRITE! 2.  905) 

905  FORMAT!/,  4X.  '4'.3X.  'SECOND  MOMENT'.  10X.  'TWO  SCANS') 

LS-0 

M2A-0  0 

M22A-0  0 

DO  120  1-1.20 

CALL  MYF ! EDGE. 1024. 1 ) 

DO  131  II-l. 1024 

131  AEDGE ! 1 1 > -EDGE ! 1 1 ) 

CALL  MYF! EDGE. 1024, 1 ) 

DO  132  11-1.1024 

132  EDGE ! 1 1 ) - ! AEDGE ! 1 1 ) ♦EDGE ! 1 1 ) ) /2.  0 
DO  137  II-LL.LH 

J-II-LL-M 

EDGE ! J ) -EDGE ! I I > 
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Attachment  2  code,  continued 


CALL  SECMO(EDGE,  til.  M2.  AVL.  AVH.  J.  K> 

IftM2>  120.  120,  129 

29  LS-LS+1 
M2A-M2A+M2 
M22A-M22A*M2**2 

WRITEO,  800 )LS,  Ml.  M2 
MR ( LS >*M2 
20  CONTINUE 

WR I TE <2. 804 >  ( J, MR< J).  J*l» LS ) 

IF(LS-2  O)  130,133.133 

33  D*SGRT  (ABS(LS*M22A-M2A**2>/ (LS*(LS— 1 ) ) > 

M2A=M2A/LS 

30  WRITE (2,  805)  M2A,  D 

AVERAGE  FOUR  EDGES  AND  COMPUTE  SECOND  MOMENT 
REWIND  1 
WRITE <2, 906) 

06  FORMAT </////.  4X.  '*',3X,  'SECOND  MOMENT'.  10X.  'FO 

LS*0 
M2A-0.  0 
M22A=0  0 
DO  220  1-1. 10 
CALL  MYF (EDGE. 1024. 1) 

DO  221  II-l. 1024 
!21  AEDGE (II) “EDGE (II) 

CALL  MYF (EDGE. 1024, 1) 

DO  222  11*1. 1024 

>22  AEDGE < II )*EDCE< II ) ♦AEDGE < II > 

CALL  MYF (EDGE, 1024, 1) 

DO  223  11*1. 1024 

?23  AEDGE (II) *EDGE (II) +AEDGE ( II ) 

CALL  MYF (EDGE.  1024,  1) 

DO  224  II-l. 1024 

224  EDGE ( 1 1 ) * ( EDGE (II) ♦AEDGE < 1 1 ) >/4  0 

DO  237  I I*LL. LH 
J-II-LL+1 

237  EDGE ( J ) *EDGE (II) 

CALL  SECMO<EDGE.  Ml.  M2,  AVL,  AVH.  J,  K) 

IF  (M2)  220.220.229 
229  LS-LS+1 

M2A-M2A+M2 
M22A*M22A«-M2**2 
MR  ( LS ) -M2 

WRITEO,  800)LS,  Ml,  M2 
220  CONTINUE 

WRITE <2. 804)  < J. MR ( J).  J-l< LS > 

IFILS-2  0)  231.232.232 

232  D-SORT ( ABS ( LS*M22A-M2A*#2 ) / <  LS* ( LS- 1 ) > ) 

M2A-M2A/LS 

231  WRITE (2, 805)  M2A,  D 

C  PERFORM  AVERAGE  OF  FORTY  SCANS  AND  RECOMPUTE  THE  M2 
REWIND  1 

WRITE <2. ») 'AVERAGE  OF  FORTY  SCANS' 

WRITE <5. *) 'AVERAGE  OF  FORTY  SCANS' 

LS*0 
M2A-0  0 


FOUR  SCANS' 
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Attachment  2  code,  continued 


M22A-0.  0 
DO  310  1*1 .  1024 
310  AEDC£(I>-0  0 

DO  321  1-1,40 

CALL  MYF (EDGE,  1024.  1 > 

DO  320  11-1, 1024 

320  AEDGE (II) -AEDCE (II) *EDCE (II) 

321  CONTINUE 

DO  331  11-1,1024 

331  EDCE (II) -AEDGE (II) /40.  0 

DO  337  II-LL.LH 
J-II-LL*1 

337  EDGE ( J ) -EDGE (II) 

CALL  SECMO ( EDGE, Ml, M2. AVL.  AVH,  J,  K> 

IF ( M2 )  341.340,340 
341  M2-ABS  ( M2 ) 

Ml  — 1000  O 
340  D-0  O 

LS-1  0 

WRITE (5. 800)  LS, Ml ,  M2 
WRITEC2,  800)  LS.M1.M2 
WRITEC2, 810)  DD.  SDD 
WR  I TE  (5,810)  DD,  SDD 

810  FORMAT ( /, 9X,  'EDGE  DENSITY  CHANGE- ', CIO.  5, SX.  'NOISE  SDEV- 

999  STOP 

END 

C  SECOND  MOMENT  SUBROUTINE 

SUBROUTINE  SECMO (EDGE, Ml.  M2,  AVL.  AVH.  J,  K) 

REAL*8  EDGE ( 1 024 ) ,  AL I NE ( 1 00 ) < X. Y.M1.M2 
REAL*8  AM 1 , AM2 
AM2—0. 0 
AM 1-0  0 

IF (AVH- AVL )  800, 800. 5 
5  DO  10  I-l.J 

10  EDGE ( I ) — ( EDGE ( I ) -AVL ) / ( AVH-AVL > 

C  CALCULATE  L(X) 

C  CALCULATE  Ml 

DO  40  1-2, J 
S-I-l 

40  AMI -AMI* (EDGE ( I )-EDGE( 1-1 ) >#S*0.  003 

WR I TE( 9, *) 'CALCULATED  Ml', AMI 
C  CALCULATE  M2  THE  FIRST  TIME 
DO  60  1-2, J 
S-I-l 

60  AM2-AM2* ( EDGE ( I )-EDGE( 1-1 > )*( (0.  003*S-AM1*0.  00i9>**2> 

GOTO  810 
800  M2— 100 

GOTO  820 
810  M2-AM2 

820  RETURN 

END 


CIO  9) 


Attachment  2  code,  continued 


SUBROUTINE  MYF(VAL, N, LF> 
INTEGER  EDGE (1 024) 

REAL* S  VAL<1024) 

CHARACTER*BO  HEADO) 

READ ( LF. 900 )  < HEAD ( I > .  1-1 .  3 ) 
FORMAT ( A80  > 

READ(LF<  901 )  < EDGE < J > .  J*1 ,  N > 
DO  10  J~1.N 

VAL (J) “EDGE ( J) / 1000.  0 
FORMAT l 2014) 

FORMAT ( IX,  2014) 

RETURN 


Read  ID 


Read  Nois 


Compute 

Density 

Variance 


Compute 

Weiner 

Spectrum 


Compute 

DQE 


Output 


Attachment  3  code 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


THE  FOR  PROGRAM  WRITTEN  BY  KONRAD  KERN 
PROGRAM  INPUTS 

1,  SECOND  MOMENT.  COMPUTED  IN  FORMER  PROGRAM 

2,  EDGE  STEP,  COMPUTED  IN  FORMER  PROGRAM 

3,  NOISE  SCAN  DATA,  105  BLOCKS  OF  64  SAMPLES 

4,  EXPOSURE  IRRADIANCE  IN  PHOTONS  PER  CM*#2 

5,  CONTRAST 

6,  SAMPLE  DATA  IDENTIFICATION 


PROGRAM  OUTPUTS: 

1.  DENSITY  VARIANCE 

2,  DQE 

3.  WEINER  SPECTRUM 

4,  INFORMATION  CONTENT 


REAL*8  D2< 105), N0ST<64) 

REAL*B  NOISE (128) 

REAL*8  I C  <  64 ) ,  MTF ( 64 ) , LS, DEN (10) , LE< 10) 
REAL*S  WS(128>. B. C.  D2BAR,  Ml- M2 
COMPLEX  NS< 128), NSC < 128) 

INTEGER  SAMPLE 
COMMON  DEN.LE 
EXTERNAL  FOUR  1 
CALL  A5SIGN ( 5,  'TI:  ') 

CALL  ASSIGN (2, 'RES  DAT') 

CALL  ASSIGN ( 3, 'NOIS  DAT') 


13  FORMAT  (F3.  2.  X,  F4.  2) 

11  FORMAT ( 20 14) 

15  FORMAT < 13) 

16  FORMAT (F5.  2) 

C  SAMPLE  IDENTIFICATION  INFORMATION 
WR I TE ( 5, * )  'SAMPLE  NUMBER  ?' 

READ <5. 15)  SAMPLE 

WRITE  <5,  #> 'N0RMAL*1.  CPA=2 ' 

READ ( 5, 15) IDS 

WRITE  (5,*)  'IMAGE  CONTRAST' 

READ (5. 12)  CON 

WRITE (5, *)  'LOG  EXPOSURE™  ? ' 

READ( 5, 12)  EXL 

WR I TE  <  5, * )  'EXPOSURE  IN  PHOTONS/AREA ' 
READ<5, 12)  QA 

C  EDGE  STEP  AND  M2  COMPUTED  IN  FORMER  PROGRAM 
WRITE (5, #)  'EDGE  STEP' 

R£AD<  5,  12>DD 
WRITE (5, *) 'TOTAL  M2' 

READ( 5, 12>M2 
AREA=731  6E-9 
FORMAT ( G1 0  5) 

WRITE (2, 802)  SAMPLE.  CON 
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Attachment  3  code,  continued 


WRITE(2. 815)  EXL 

815  FORMAT (/, 5X. 'LOG  EXPOSURE* F8  3 ) 

802  FORMAT ( /.  5X.  'SAMPLE  NUMBER-'.  13.  10X.  'CONTRAST- '.  F5  3) 
WR  ITE (2. 803 )  QA.  IDS 

WRITE (2. 850  > M2 

850  FORMAT  (  5X.  'M2  USED*',  CIO.  5) 

803  FORMATt /. 5X. 'EXPOSURE*  ', G1 1  5, 5X, 'EXPOSURE  TYPE*  '.1 
C  BEGIN  COMPUTING  THE  NOISE  INFORMATION 

DO  400  1-1.  128 
WS(I )*0  0 
400  D2( I )*0. 0 

D2BAR-0  0 
DO  540  J=l,  105 
CALL  MYF(NOST,  64,  3) 

B-0  00 
C-0  00 

C  COMPUTE  DENSITY  VARIANCE 
DO  550  1-1. 64 

NO I SE ( I ) =NOST  < I ) 

B-B+NOST  < I ) 

C-C+NOST ( I )*NOST( I ) 

550  CONTINUE 

D2(J)*ABS(64  0*C-B*B ) / ( 64.  0*63.  0> 

D2BAR=D2BAR*D2< J) 

WRITE (5. 901 ) J. D2( J),  D2BAR 

901  FORMAT (IX#  15.  2X,  'D2*  IX.  Oil.  5.  3X.  'D2BAR*  '»  IX, Gil  5> 

C  REMOVE  DC  LEVEL 
A=0  O 

DO  570  1-1. 64 
570  A*A+NOISE(I > 

A=A/64. 0 
DO  580  1-1. 64 

580  NO I SE ( I  ) -NO I SE ( 1 ) —A 

C  APPLY  BARTLETT  WINDOW 
DO  590  1*1.32 

590  NS ( I ) — CMPLX ( NOISE ( I >  *  C 1 . 0-1. 0/32  0*1 ) ) 

DO  591  1-1,32 

591  NS( 129-1 ) —CMPLX ( NOISE ( I *32 )* ( 1 .  0-1.  0/32  0*1 ) ) 

DO  600  1-33.  96 

600  NS ( I )*(0  0, 0. 0) 

C  TRANSFORM  TO  FREQUENCY 

CALL  F0UR1 (NS. 128, 1 ) 

DO  610. 1*1, 128 

NSC ( I ) — CON JG  <  NS ( I ) ) 

610  NOISE < I ) -(CABS (NS ( I ) *NSC ( I >  > >**2 

C  PERFORM  ENSEMBLE  AVERAGE  OF  WEINER  SPECTRUM 
DO  630  1-1. 128 

630  WS ( I ) “WS  < I ) +N0 1 5E ( I ) 

540  CONTINUE 

C  AVERAGE  WEINER  SPECTRUM  AND  PERFORM  FFT  SCALING 
DO  640  1*1. 128 

640  WS( I >-WS( I >/105  0/128  0 

D2BAR=D2BAR/ 105  0 
SDD-SQRT (D2BAR ) 

AL2-AL0G ( 2  0) 
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Attachment  3  code,  continued 


C  COMPUTE  FIRST  MOMENT  OF  WEINER  SPECTRUM 
M1*0  0 

DO  630  1=2.64 
S=I-1 

650  M 1 =M 1 ♦ ( WS < I ) -WS ( I - 1 )  )*S*2.  604 

WRITE(2. 606)  Ml 

606  FORMAT(/,  5X.  'WEINER  SPECTRUM  Ml  =  '.Cll. 3) 

WRITE <2, 903) 

903  FORMAT!/,  10X,  'WEINER  SPECTRUM'./) 

WRITE (2. 902)  (WS ( I ) . 1=1 . 64 ) 

902  FORMAT ( IX,  Cll.  5.  SX.  CM.  5.  5X,  Cll.  5) 

C  COMPUTE  INFORMATION  CONTENT 
SUM=0. 0 
B=1  666**2 
A=DD**2/4. 0 
DO  660  1=1.64 

F*FLOAT( 1-1 )*2.  604 
MTF ( I )*1— 2. 0*3.  1416**2*F**2*M2 
IF  < MTF <  I  )  )  636.656,658 
656  MTF  <  I  )  =0.  O 

658  P«A/ (B«-F**2> 

0*1. 0+P*(MTF< I >/WS< I ) ) 

660  IC< I )=AL0C(0)/AL2*F 

C  INTEGRATION  BY  SIMPSON'S  RULE 
DO  670  1=2.  64 

670  SUM* SUM* ( IC ( I ) *IC ( I—l ) ) /2.  0*2  604 

SUM=SUM*3  1416 
C  OUTPUT  CALCULATED  VALUES 
WR1TE12, 610)  DD. SDD 

BIO  FORMAT ( /, 5X>  'EDGE  DENSITY  CHANGE- '. CIO  3. 5X. 'NOISE  SDEV='.G10  5) 

WRITE <2. 607)  SUM 

B07  FORMAT (/, 3X.  ' INFORMATION  CAPACITY  =  '.Cll,  5) 

DQE=0  1 B86*C0N*C0N/ <  GA*AREA*D2BAR ) 

950  FORMAT ( C16  5) 

WRITE <2. 808)  AREA, DEBAR,  DOE 

808  FORMAT (/,  2X,  'AREA-',  Cll,  3,  3X,  'D2=',  Cll.  3.  3X.  'DOE-',  Cll.  3.  ///) 

999  STOP 

END 

SUBROUTINE  MYF ( VAL.  N, LF ) 

INTEGER  I EDGE ( 1024 ) 

REAL*B  VAL (1024) 

CHARACTER*80  HEAD ( 3  > 

READ (LF, 900)  < HEAD ( I ) .  1  =  1 . 3 ) 

900  FORMAT (A80) 

READ(LF. 901 )  < 1 EDGE < J ) . J= 1 . N ) 

DO  10  U*1.N 

VAL  <  J )  =  I EDGE  <  J ) 

10  VAL ( J ) » VAL  <  J ) / 1 000 .  0 

901  FORMAT (2014) 

902  FORMAT ( IX, 2014) 

RETURN 

END 
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